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detonated  underwater.  The  plates  were  instrumented  to 
measure  transient  strains  and  free  field  pressures.  Test 
results  were  compared  to  preshot  calculations  done  using  the 
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INTRODUCTION 


I . 

A.  BACKGROUND  AND  HISTORY 

Sergei  Gorshkov,  Admiral  of  the  Soviet  Fleet,  has 
observed  that  the  end  of  World  War  Two  marked,  "... 
the  start  of  a  military-technical  revolution  which  in  scope 
and  depth  transcended  all  the  reforms  and  transmissions 
which  had  previously  occurred  in  the  armies  and  fleets  of 
the  world."  [Ref.  1] 

Much  has  been  made  of  the  tactical  and  geo-political 
implications  of  that  revolution.  Less  popularly  publicized, 
but  equally  radical  are  the  changes  that  the  developments  in 
te-chnology  have  spawned  in  the  processes  of  ship  design  and 
construction . 

Throughout  most  of  history,  war  at  sea  has  been  a 
•'*  '  of  trying  to  set  one's  opponent  afire,  board  him,  or 

aurl  chu. ks  of  stone  or  metal  at  him  across  relatively  short 
distances.  All  of  these  methods  had  the  distinguishing 
characteristic  of  being  essentially  disabling  tactics, 
unless  a  ship  was  burned  or  perforated  to  the  point  of 
sinking  (and  they  often  were  not),  damage  was  generally 
largely  superficial  and  the  ship  could  be  refitted  to  see 
service  again.  Witness  the  many  vessels  that  have  been 


fought  on  both  sides  of  a  conflict  and  the  tradition  of 


incident  wave  is  impinging  upon  a  medium  of  greater  density 
than  the  fluid.  The  reflected  wave  will  therefore  be 
positive.  The  actual  magnitude  of  this  wave  depends  upon 
the  physical  characteristics  of  the  bottom  material.  The 
effects  of  this  bottom  reflected  wave  are  also  additive 
and  will  generally  reach  the  target  at  some  time  after 
the  initial  shock  front  has  passed;  although  the  actual 
time  of  arrival  of  this  secondary  front  will  once  again 
be  dependent  upon  the  bottom  conditions  and  the  speed  at 
which  the  wave  travels  through  the  bottom  medium.  Addi¬ 
tionally,  in  shoal  waters,  that  same  bottom  reflected  wave 
will  reach  the  surface;  creating  its  own  surface  cutoff. 

The  net  effect  upon  the  target  is  thus  one  of  an  inci¬ 
dent  shock  wave;  which,  if  the  charge  and  target  are  suffi¬ 
ciently  near  the  surface,  will  be  cutoff  by  the  negative 
surface  wave.  This,  in  turn,  is  followed  by  the  bottom 
reflected  wave  which  also  experiences  an  exponential  decay, 
and  which  may  also  be  cutoff  by  its  own  generated  surface 
wave  ( Fig .  2.3). 

Analytically,  the  easiest  way  to  model  the  effects  of 
incident,  surface,  and  bottom  waves  upon  the  target  is  as 
has  been  represented  in  Figure  2.4,  with  three  separate 
images . 

The  incident  and  surface  images  are  relatively  easy  to 
depict.  The  incident  image,  ,  of  course  represents  the 
charge  itself  and  occupies  the  same  location  in  the  spatial 
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referred  .to  as  surface  cutoff,  and  the  negative  pressure  as 
bulk  cavitation. 


Figure  2.2  Exponential  Approximation  with  Surface  Cutoff 


The  direct  effect  of  bulk  cavitation  is,  by  itself,  not 
particularly  harmful.  The  structure  simply  experiences 
essentially  no  pressure  while  it  is  in  the  cavitated  region 
When,  however,  the  static  head  of  water  above  the  region 
overcomes  the  cavitation  forces,  the  region  closes  suddenly 
to  generate  a  destructive  reloading  upon  the  structure. 

To  further  complicate  matters,  if  the  shot  occurs  in 
shoal  waters  or  at  depths  near  the  bottom  a  second  boundary 
is  also  present.  Unlike  the  surface,  in  this  case  the 


compliance  with  the  laws  of  continuity.  This  implies  that 
a  negat ive  shock  wave  of  magnitude  equal  to  the  original 
incident  wave  must  be  propagated  downward  through  the 
medium  to  satisfy  equilibrium  conditions.  This  negative 
wave  will  travel  through  the  water  at  an  angle  to  the 
surface  that  is  equal  to  that  of  the  incident  wave;  much 
as  the  negative  reflection  image  from  a  mirror  continues  at 
an  angle  equal  to  the  angle  of  incidence. 

Remembering  that  these  waves  are  radiating  spherically 
outward  from  the  point  of  inception,  it  is  clear  that  both 
the  incident  wave  and  the  generated  surface  wave  will  even¬ 
tually  reach  the  ship,  submarine,  or  other  target  object. 

It  is  likewise  clear  that,  with  its  longer  path,  the  surface 
wave  must  reach  the  target  at  some  time  later  than  the  inci¬ 
dent;  which  travels  only  the  shortest,  or  standoff,  distance 
between  charge  and  target.  The  net  effect  is  additive. 

When  the  surface  wave  reaches  the  target,  the  initial  shock 
front  will  have  passed  and  the  incident  wa\ e  will  be  at  some 
point  in  its  exponential  decay.  The  negative  surface  wave 
will  then  "chop  off"  the  tail  of  the  incident  wave.  If  the 
spatial  geometry  of  the  shot  is  such  that  the  magnitude  of 
the  surface  wave  is  greater  than  the  remaining  magnitude  of 
the  incident  wave,  a  region  of  negative  pressure  will  be 
created  (Fig.  2.2).  If  this  negative  is  less  than 
the  pressure  required  to  keep  the  fluid  in  a  liquid  state, 
the  fluid  will  "flash"  into  vapor.  This  chopping  effect  is 
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where  PQ  is  in  pounds  per  square  inch  and  f  in  milli¬ 
seconds.  ,  Kg ,  ,  Ag  are  the  explosive  dependent  values 

(Table  I).  W  and  R  are  the  charge  weight  (in  pounds)  and 
standoff  distance  (in  feet),  respectively. 

Table  I 

Explosive  Constants  I 


Explosive 

TNT 

HBX-1 

PENTOLITE 

K1 

22505 

22347.6 

24589 

A1 

1.18 

1.144 

1.194 

K2 

.058 

.056 

.052 

A0 

-.185 

-.247 

-.257 

B.  SURFACE  CUTOFF  AND  BOTTOM  REFLECTION 

Shock  waves  propagate  radially  from  the  charge.  Thus, 
for  ‘a  point  source  in  an  infinite  medium,  the  wave  could  be 
expected  to  travel  as  an  ever  expanding  sphere  until  damped 
by  the  fluid.  For  small  to  moderately  sized  charges  deto¬ 
nated  well  below  the  surface  in  great  depths  of  water,  this 
is  a  satisfactory  analogy.  For  explosions  nearer  the 
surface  or  in  shoal  water,  however,  there  are  secondary 
effects  for  which  there  must  be  an  accounting. 

When  a  shock  wave  reaches  the  surface  of  the  water,  the 
adjoining  air  is  of  sufficiently  lesser  density  to  be  inca¬ 
pable  of  supporting  the  shock  wave.  The  result  is  that  the 
effective  pressure  at  the  interface  must  be  zero  in 


functions  of  the  type  of  explosive,  the  charge  weight,  and 
the  standoff  distance.^"  The  equations  for  the  pressure  and 
decay  constant  were  originally  determined  by  Robert  H. 

Cole  and  were  presented  in  his  work,  Underwater  Explosions 
[Ref.  5].  These  have  since  been  simplified  to  use  a  series 
of  empirically  derived  coefficients  developed  by  R.  S. 

Price  [Ref .  7] . 

The  pressure  profile  can  be  expressed  by  the  equation 
t1  -  t 

P(t)  =  PQexp  - q -  (eqn  2.1) 


where  Pq  is  the  initial  (greatest)  pressure  of  the  shock 
wave,  0  is  the  decay  constant  describing  the  exponential 
decay,  and  t  -  t^  is  the  time  elapsed  from  the  arrival  of 
the  shock  front. 

The  initial  pressure,  Pq ,  and  the  decay  constant, 

0,  are  dependent  largely  upon  the  type  of  explosive  and 
the  weight  of  the  charge.  These  may  be  expressed 


P 

o 


K. 


|V/3y2 


and , 
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V1,3(^)42 


( eqn  2.2) 


( eqn  2.3) 


The  standoff 
the  charge  to  the 


distance  is  defined  as  the  distance  from 
nearest  point  on  the  target . 


0 


•  .*•  *'•  .* 
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therefore  may  be  apprehended  as  an  essentially  acoustic 
phenomenon.  As  the  water  begins  to  be  displaced,  the  gas 
bubble  created  by  the  explosion  expands;  thereby  decreasing 
its  internal  pressure  and  correspondingly  the  pressure  on 
the  constraining  fluid.  Viewed  on  a  pressure-time  curve  the 
initial  rise  in  pressure  is  so  nearly  vertical  as  to  appear 
discontinuous,  while  the  pressure  relief  is  characterized  by 
an  approximately  exponential  decay  (Fig.  2.1)  [Ref.  5:  pp . 
4-7]. 


II. 


UNDERWATER  EXPLOSION  THEORY 


A.  THE  PRESSURE  PULSE 

Any  conversion  of  matter,  whether  chemical  or  nuclear, 
which  results  in  the  very  rapid  production  of  large  quanti¬ 
ties  of  gas  at  very  high  temperature  can  be  categorized 
within  the  broad  species  "explosion.”  Typically,  for  mili¬ 
tary  high  explosives  such  as  TNT  the  pressure  and  tempera¬ 
ture  within  the  gases  produced  by  an  explosion  are  on  the 
order  of  50,000  atmospheres  and  3000  degrees  Celsius 
[Ref.  5:  p.  3].  For  nuclear  explosions  the  initial  tempera¬ 
ture  is  on  the  order  of  a  million  degrees  Celsius  and  the 
pressure,  which  is  dependent  upon  the  yield  is  similarly 
orders  of  magnitude  higher  than  that  for  a  conventional 
explosive  [Ref.  6]. 

Clearly,  temperatures  and  pressures  of  the  levels 
described  must  dissipate  high  levels  of  energy  through  any 
surrounding  medium.  In  the  case  of  water,  the  initial  mani¬ 
festation  of  this  disturbance  is  an  intense  compression 
created  pressure  wave  propagating  radially  outward  from  the 
charge;  followed  nearly  instantaneously  by  a  displacement  of 
the  fluid.  In  the  immediate  vicinity  of  the  charge,  the 
velocity  of  propagation  of  the  pressure  wave  is  several 
times  the  speed  of  sound  (about  5000  feet  per  second  in 
water)  but  approaches  the  value  almost  immediately  and 
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retained.  Conclusions  regarding  the  ability  of  USA-STAGS  t 
predict  the  responses  of  a  simple  model  to  underwater  explo 
sions  will  be  drawn;  and,  perhaps  of  even  greater  impor¬ 
tance,  guidelines  will  be  established  for  using  the 
USA-STAGS  code  and  conducting  underwater  explosion  testing 


for  future  studies. 


conditions  that  could  be  expected  in  a  ship's  plating;  and 
were  subjected  to  shock  loadings  that  produced  deformations 
well  into  the  plastic  regime  in  order  to  ensure  stiffener 
tripping.  Existing  information  concerning  underwater  explo 
sion  phenomena  was  employed  in  an  attempt  to  reduce  experi¬ 
mental  uncertainties. 

The  first  study  in  this  series  was  that  conducted  by  LT 
Thomas  Rentz,  USN  in  fulfillment  of  thesis  requirements  at 
the  Naval  Postgraduate  School  [Ref.  3]  and  [Ref.  4].  In 
that  study  he  developed  the  following  experimental  strategy 
to  study  the  EPSA  code: 

An  attempt  to  validate  the  code  was  conducted  in  two 
phases:  Using  pressure  approximations  based  upon  empiri¬ 

cally  determined  formulae,  the  code  was  first  used  as  a 
pre-test  predictor  of  test  results.  Once  the  underwater 
explosion  testing  had  been  conducted,  pressure  histories 
derived  from  that  testing  were  input  into  the  code.  The 
plate  strains  predicted  by  the  code  in  both  its  preshot  and 
post-shot  capacities  were  then  compared  to  actual  strain 
histories  from  testing.  Plate  deformations  and  boundary 
stresses  predicted  by  the  code  were  compared  to  observed 
Dhvsical  responses. 

In  this  study,  the  code  to  be  tested  is  USA-STAG3.  In 
an  effort  to  maintain  continuity  with  the  earlier  study  in 
the  series,  the  same  experimental  strategy  will  be  followed 
The  test  geometry  and  flat  plate  model  will  also  be 


What  has  been  needed  is  a  method  of  analyzing  the  design 
itself  in  order  to  achieve  a  numerical  approximation  of  the 
end  product's  response  to  shock. 

This  need  has  coalesced  into  the  development  of  two 
computer  codes;  EPSA  (Elasto  Plastic  Shell  Analysis), 
produced  by  Weidlinger  Associates,  and  USA-STAGS  (Underwater 
Shock  Analysis — Structural  Analysis  of  General  Shells) 
produced  by  Lockheed  Missiles  and  Spacecraft,  Inc. 

B.  OBJECTIVES 

This  study  is  the  second  in  an  ongoing  series  sponsored 
by  the  Defense  Nuclear  Agency  (DNA)  into  the  effects  of 
underwater  explosions  upon  ship  hulls.  The  intent  of  these 
studies  has  been  to  test  the  applicability  of  the  two  previ¬ 
ously  mentioned  structural  response,  underwater  shock  anal¬ 
ysis  computer  codes  to  naval  design  and  analysis  requirements 
Under  the  umbrella  of  general  applicability,  focus  has 
been  placed  upon  two  defined  goals: 

1.  Use  of  the  codes  to  gain  insights  into  large  deflec¬ 
tion  elastic-plastic  responses  of  a  submerged  struc¬ 
ture  subjected  to  transient  acoustic  shock  loading; 
with  especial  emphasis  upon  stiffener  tripping 
phenomena . 

2.  Performance  of  underwater  explosion  testing  of  the 
structure  to  validate  the  results  predicted  by  the 
code  . 

To  achieve  these  goals,  flat  plate  models  geometrically 
similar  to  a  ship's  stiffened  hull  structure  have  been 
studied.  These  plates  were  air  backed  to  simulate  general 


withstand  specified  levels  of  shock  loading.  These  re¬ 
quirements  included  (and  continue  to  include)  actual  shock 
testing  of  the  first  ship  of  a  class  and  many  crucial  systems 
on  the  component  level. 

During  the  immediate  post-World  War  Two  period,  many  of 
these  requirements  were  largely  tentative  and  experimental. 

By  the  early  1960's,  however,  sufficient  experience  had  been 
gained  to  allow  design  requirements  to  be  set  forth 
[Ref.  2].  These  specif icat ions  established  specific  guide¬ 
lines  for  contractors  to  meet  based  upon  type  of  equipment, 
type  of  vessel  in  which  the  equipment  was  to  be  used,  and 
location  and  method  of  mounting  within  the  vessel.  As  with 
the  ships  themselves,  actual  shock  testing  of  the  item  was 
the  primary  means  of  establishing  compliance;  as  it  remains 
today . 

Performing  a  shock  test  upon  production  end  products, 
particularly  when'  those  items  are  as  large  as  an  aircraft 
carrier  can  have  some  obvious  drawbacks:  it  is  costly,  it 
is  dangerous,  and  the  results  can  be  embarrassing. 

Moreover,  there  can  be  longterm  structural  damage  which 
is  neither  immediately  apparent  nor  correctable,  but  which 
does  not  by  itself  constitute  non-compliance  of  contract 
terms.  Moreover,  some  items  simply  do  not  lend  themselves 
to  shock  testing. 


generation  of  intense  pressure  waves  in  the  sea  that  could 
break  a  vessel's  back  or  rupture  hull  fittings  whether  the 
vessel's  skin  was  punctured  or  not.  ''imply  increasing  the 
thickness  of  plating  in  discrete  areas  of  vulnerability  was 
no  longer  sufficient.  Design  for  strength  and  resiliency 
under  unequal  pressure  loadings  was  required.  Moreover, 
this  requirement  applied  to  attacker  as  well  as  attacked. 
Although  the  attacking  surface  vessel  had  the  advantage  of 
being  able  to  run  from  the.  immediate  scene  of  the  explosion 
and  the  effects  were  not  compounded  by  as  great  a  pre¬ 
existing  static  head  of  water,  these  weapons  were  not  direc¬ 
tional  and  could  be  expected  to  damage  all  within  their 
sphere . 

The  issue  of  a  vessel's  ability  to  withstand  an  under¬ 
water  explosion,  be  it  submarine  or  surface  ship,  became 
more  pressing  as  the  explosive  devices  became  increasingly 
powerful  and  sophisticated.  Not  only  did  depth  charges 
become  lethal  over  a  wider  area,  but  the  advancing  tech¬ 
nology  permitted  mines,  too,  to  grow  from  essentially 
contact  weapons  to  devices  that  could  endanger  all  types  of 
shipping  from  considerable  depths.  The  issue  has  become 
particularly  crucial  with  the  development  of  nuclear  weapons 
which  are  so  powerful  that  no  device  need  contact  a  ship  or 
group  of  ships  to  cause  rampant  destruction. 

The  response  of  the  U.S.  Navy  was  to  require  that  ships 
and  their  critical  components  be  designed  and  built  to 
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expecting  naval  personnel  to  supplement  their  pay  with  prize 
money . 

Two  events  occurred  at  or  around  the  beginning  of  the 
twentieth  century  to  change  all  of  that:  the  invention  of 
high  explosives  and  the  advent  of  the  submarine. 

Initially,  the  submarine  was  primarily  a  commerce 
raider.  It  was  too  small  to  carry  a  prize  crew  and  too 
lightly  armed  to  effectively  combat  a  conventional  warship. 
During  World  War  One,  most  battles  between  submarines  and 
conventional  ships  were  fought  on  the  surface.  Submarines 
would  use  their  stealth  to  sneak-up  upon  unsuspecting  mer¬ 
chantmen,  surface  and  then  attempt  to  sink  them  with  their 
deck  cannons.  Torpedoes  were  largely  unreliable  and  the 
submarines  themselves  were  at  risk  if  fast  armed  vessels 
were  present.  Convoying  the  merchants  with  small,  agile 
men-of-war  was,  briefly,  an  effective  means  of  countering 
the  submarine  threat.  However,  as  torpedoes  became  more 
effective  so,  too,  did  the  submarine.  For  the  convoy  escort 
vessels  it  was  no  longer  sufficient  to  have  to  wait  for  a 
submarine  to  surface  in  order  to  attack  it.  An  effective 
standoff  weapon  was  needed. 

That  standoff  weapon  was  the  depth  charge.  First  intro¬ 
duced  in  World  War  One,  the  depth  charge  came  of  age  in  the 
Second  World  War;  and  with  it  a  new  era  in  ship  design  was 
ushered  into  existence.  For  the  principle  of  this  new 
weapon  was  not  one  of  punching  a  hole  in  a  ship's  hull  and 
then  exploding  internally;  rather,  high  explosives  permitted 
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is  of  the  same  magnitude 


geometry.  The  surface  image,  il  , 
as  the  incident  image  but  of  an  opposite  sign.  It  lies  an  a 
line  drawn  vertically  through  the  charge  at  a  distance  above 
the  surface  equal  to  the  depth  of  the  charge  below. 

The  surface  cutoff  pressure  may  therefore  be  computed 
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where  R  is  now  the  distance  from  the 
s 

the  target.  The  time  of  cutoff  after 
initial  shock  front  i‘s  thus  the  time 
to  cover  the  greater  distance,  or 


( eqn  2.4) 

surface  image  to 
the  passage  of  the 
difference  for  the  wave 


t 


c 


( eqn  2.5) 


where  c  is  the  speed  of  sound  in  the  fluid.  .The  cutoff 
pressure  will  then  decay  in  the  same  exponential  manner  as 
the  incident  pressure. 

As  may  be  conceived  from  the  preceding  discus-ion  of  the 
bottom  reflection;  the  location  of  the  bottom  image  is  not 
so  straightforward.  It  may  be  viewed  as  initiating  along 
the  same  vertical  line  as  do  the  incident  and  surface 
images,  and  at  a  depth  below  the  fluid-bottom  interface 
equal  to  the  distance  of  the  charge  above.  Here,  however, 
the  similarities  to  the  surface  image  end.  The  bottom  image 
cannot  be  considered  a  point  charge  because  of  the  tendency 
of  the  pressure  wave  to  "smear"  across  the  bottom.  The 
magnitude  and  speed  of  the  wave,  too,  are  dependent  upon 
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bottom  characteristics.  A  muddy  bottom,  for  instance,  will 
tend  to  reflect  a  diffuse,  low  magnitude  wave  travelling  at 
or  about  the  speed  of  sound  in  water.  A  rock  bottom  on  the 
other  hand  will  produce  sharp,  high  magnitude  pulses  which 
will  tend  to  radiate  over  a  distance  due  to  the  higher  speed 
of  travel  of  the  shock  wave  in  the  bottom  medium  than  in  the 
fluid.  Due  to  this  dependence  upon  bottom  conditions, 
determination  of  bottom  reflection  effects  do  not  lend  them¬ 
selves  to  general  predictive  studies  and  are  often  disre¬ 
garded,  although  they  can  be  significant. 


C.  BUBBLE  EFFECT 

The  foregoing  discussion  of  the  development  of  the  pres¬ 
sure  wave  and  its  impact  upon  the  target  assumes  that  a 
single,  clean  wavefront  is  produced  by  the  explosion.  In 
point  of  fact,  such  is  rarely  the  case.  When  the  charge 
explodes,  it  first  creates  a  high  intensity  shock,  front 
which  is  the  initial  incident  shock  wave.  It  also  produces 
a  bubble  of  very  hot  gasses  at  very  high  pressure. 

The  bubble  expands  until  it  reaches  pressure  equilibrium 
with  the  surrounding  water.  When  it  does  so,  a  portion  of 
the  energy  within  the  bubble  is  released  as  a  new  shock 
front.  This  release  of  energy  then  permits  the  momentum  of 
the  water  displaced  to  assert  itself  and  the  bubble  is 
collapsed,  compressing  the  gases  within.  The  displaced 
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pre-bubble  condition. 


As  the  water  rushes  in,  it  compresses 


The  time  and  radius  of  the  bubble  at  the  first  pulse 
has,  again,  been  determined  empirically,  the  maximum  radius 


max  3^(D  +  33)1/2 


( eqn  2.6) 


at  time, 


t  =  K 


H(D  +  33 )°  /S 


( eqn  2.7) 


where  A  is  expressed  in  feet  and  t  in  seconds.  W  is 
max 

once  more  the  weight  of  the  charge,  and  are  explosive 
dependent  constants  (Table  II),  and  D  is  the  depth  of  the 
charge  in  feet . 


Table  II 


Explosive  Constants  II 


Explosive 


HBX-1 


PENTOLITE 


14.14 


12.67 


12.88 


4.761 


4.268 


4.339 


A  second  and  distinct  phenomenon  is  bubble  migration;  s 
called  because  of  the  tendency  of  the  bubble  to  migrate 
toward  any  nearby  mass.  When  the  bubble  contacts  the  mass, 
it  collapses  and  begins  a  series  of  rapid  pulsations.  The 
net  effect  is  similar  to  that  of  the  already  described 
bubble  pulse;  a  series  of  reloadings  upon  the  structure. 


The  rapidity  of  the  pulsations  of  the  collapsed  bubble  and 


their  increased  localization  tend,  however,  to  magnify  their 
severity . 

D.  APPLICATIONS  OF  THEORY  TO  EXPERIMENT  DESIGN 

Developing  a  basic  understanding  of  the  primary 
phenomena  involved  in  underwater  explosion  theory  is 
crucial  to  developing  a  successful  underwater  test  program 
for  many  reasons.  Two  are  of  particular  interest  here, 
however : 

First,  and  most  obviously,  equations  2.1  through  2.3 
allow  initial  estimates  of  the  peak  pressure  and  pressure 
history  to  which  a  test  structure  is  to  be  exposed.  These 
are  then  input  to  the  computer  code  so  that  preshot  approxi¬ 
mations  of  the  structural  responses  may  be  made;  thereby 
creating  guidelines  for  calibrating  gages,  determining 
charge  weights,  etc. 

Second,  the  design  of  the  test  itself  is  a  direct  conse¬ 
quence  of  the  phenomena  being  observed.  As  is  readily 
apparent  from  the  earlier  sections  of  this  chapter,  the 
loadings  imposed  upon  a  structure  by  an  underwater  explosion 
can  be  extremely  complex.  The  incident  wave  may  be  a 
simple  exponential  decay  dependent  only  upon  charge  and 
standoff,  but  the  many  secondary  effects  are  subject  to 
variables  which  cannot  be  readily  controlled  experimentally. 

The  contributions  of  the  secondary  effects  are  signifi¬ 
cant  and  must  eventually  be  studied.  However,  during  these 
early  phases  of  investigation  into  underwater  explosion 


effects  a  determined  effort  has  been  made  to  simplify  the 
problem  as  much  as  possible. 

It  has  therefore  been  required  that: 

1.  The  model  to  be  studied  be  as  basic  as  is  consis¬ 
tent  with  stated  objectives. 

2.  Study  be  limited  only  to  the  incident  shock  wave. 

These  two  requirements  are  interrelated  and  are  tied 

intrinsically  to  the  nature  of  underwater  explosion 
phenomena.  As  has  been  seen  earlier,  the  incident  shock  may 
be  approximated  as  a  plane  wave  emanating  from  the  charge 
and  striking  the  panel.  In  the  absence  of  secondary 
effects,  the  responses  being  studied  here  may  thus  be 
limited  to  those  of  the  test  panel  alone.  If,  on  the  other 
hand,  secondary  effects  are  included,  more  than  a  single 
point  source  is  involved  and  the  interactions  of  any  support 
or  backing  structures  with  the  fluid  medium  must  also  be 
considered. 

Fortunately,  for  the  purposes  of  this  study,  surface 
cutoff  and  bottom  reflection  are  later  time  responses  which 
can  be  neglected  if  study  is  confined  to  the  first  few 
milliseconds  after  detonation.  Bubble  and  cavitation 
effects  can  be  reduced  greatly  or  eliminated  through  careful 
attention  to  the  test  geometry. 

If  an  explosion  occurs  sufficiently  near  the  surface 
that  the  initial  bubble  radius  is  greater  than  the  charge 
depth,  the  hot  gases  forming  the  bubble  will  vent  into  the 
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Bulk  cavitation  and  bubble  effects  can  therefore  be 
reduced  significantly  or  eliminated  entirely  if  the  test 
geometry  is  such  that  the  charge  is  located  directly  above 
the  plate  at  a  depth  to  cause  bubble  venting. 

In  designing  the  test  platform  [Ref.  3:  pp .  19-28]  a 
series  of  iterations  were  run  using  the  EPSA  code  to  deter¬ 
mine  a  charge  weight  and  standoff  distance  that  would 

produce  maximum  test  panel  deflections  of  approximately  four 

2 

plate  thicknesses.  From  these  runs  it  was  determined  that 
the  ideal  charge  weight  and  standoff  w,-re  8  pounds  of  TNT  at 
9  feet. 

The  maximum  bubble  radius  is  expressed  as  a  function  of 
charge  weight  and  depth  (Eqn.  2.6).  Since  it  has  been  found 
that  bubble  venting  occurs  when  the  ratio  of  charge  depth  to 

3 

maximum  radius  is  less  than  0.75,  choosing  a  charge  depth 
is  a  fairly  straightforward  procedure  once  a  determination 
of  the  desired  ratio  has  been  made. 

2 

This  value  was  recommended  by  Weidlinger  Associates  as 
one  which  experience  indicated  would  ensure  stiffener 
tripping. 

3 

Any  ratio  below  0.75  will  cause  venting.  Venting  by 
its  very  nature  will  duct  a  portion  of  the  .explosion  energy 
into  the  atmosphere,  thereby  decreasing  the  maximum  pressure 
of  the  shock  wave.  As  a  rule  of  thumb,  this  decrease  in 
usually  considered  negligible  until  the  ratio  drops  below 
0.50. 


By  selecting  a  depth  to  radius  ratio  of  0.50  it  was 
readily  determined  that  a  charge  depth  of  four  feet  was 


appropriate.  Applying  this  to  Eqn .  2.6  yields 


A 


max 


(12.67) 


8 


1/3 


(4  +  33) 


=  7.60  ft 


4  ft 
7.60  ft 


0.526 


which  assures  venting. 

From  the  foregoing  it  can  be  seen  that  a  test  configura¬ 
tion  that  will  have  the  desired  consequences  of  reducing 
bubble  and  bulk  cavitation  effects  can  be  designed.  For  an 
eight  pound  TNT  charge,  this  platform  will  have  the  charge 
centered  vertically  at  a  standoff  distance  of  nine  feet 
above  the  plate.  The  charge  itself  will  be  located  at  a 
depth  of  four  feet  below  that  surface.  The  actual  design  of 
the  test  platform  will  be  discussed  more  fully  in  Chapter  V. 


III.  THE  FLAT  PLATE  MODEL 


As  stated  in  the  introduction,  the  objective  of  this 
study  has  been  to  test  the  applicability  of  a  computer  code, 
USA-STAGS,  to  predicting  the  effects  of  an  underwater  explo¬ 
sion  upon  ships'  hulls.  To  do  this,  it  is  clearly  necessary 

that  some  basis  for  modelling  simple  sections  of  a  typical 
4 

ship's  grillage  be  established. 

The  scaling  laws  used  here  and  in  the  earlier  work  of 
this  series  [Ref.  3]  were  those  developed  by  Dr.  Raymond  P. 
Daddazio  of  Weidlinger  Associates,  Inc. 

In  his  work  [Ref.  8],  Dr.  Daddazio  has  established  two 
dimensionless  parameters,  3  and  X,  which  he  refers  to  as 
the  "plate  slenderness  ratio,”  and  the  "longitudinal  stif¬ 
fener  slenderness  ratio,”  respectively.  These  are  defined 

,  _  b  °_y_  ( eqn  3 . 1 ) 
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and 


(eqn  3.2) 
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The  ship's  hull 
ated  stiffeners. 


or  superst ruct ural 


plating  and  associ- 


where 


a  =  material  vield  strength 

y 

E  =  Young's  modulus 

a  =  spacing  between  transverse  frames 

b  =  spacing  between  longitudinal  stiffeners 

t  =  plate  thickness 

k  =  radius  of  gyration  of  the  longitudinal  stiffener 

acting  with  an  assumed  effective  width  of  plating 

b  where, 
e 


Through  sampling  a  variety  of  ship  types,  Dr.  Daddazio 
found  that  typical  values  for  ships’  grillage  fall  into  the 
range : 

1-.0  <  3  <  4.5 
and , 

0.15  <  X  <0.9 

The  intent  of  this  series  of  studies  has  been  to  use  one 
basic  model  but  with  a  variety  of  stiffener  types  so  that 
the  responses  of  different  stiffener  geometries  could  be 
studied. 

The  model  established  in  Ref.  3  was  a  3/16  inch  thick 
test  panel  18  inches  long  by  12  inches  wide  milled  into  the 
center  of  a  two  inch  thick  6061-T6  aluminum  blank  measuring 


33  inches  by  27  inches.  Two  stiffeners  with  rectangular 
cross  sections,  3/16  inch  thick  by  one  inch  deep  were  mounted 
transversely  across  the  test  panel  in  such  a  manner  as  to  be 
symmetric  to  both  the  transverse  and  longitudinal  axes  (Fig. 

3.1) .  The  entire  assembly  was  machined  as  a  single  unit  to 
avoid  the  inconsistencies  of  welding.  A  radius  was  cut  at 
all  corners  and  the  whole  polished  to  avoid  stress  concen¬ 
tration  points. 

To  apply  Equation  3.1  and  Equation  3.2  the  two  milled 
stiffeners  were  considered  to  be  longitudinal  stiffening 
members  and  the  longer  boundaries  of  the  panel  to  be 
transverse  frame  members.  Thus,  the  dimensions  of  the  plate 
are,  a  =  12  inches,  b  =  6  inches,  and  t  =  3/16  inches.  Using 
the  standard  reference  values  for  6061-T6  aluminum,  it  was 
found  that,  for  this  original  plate:  [Ref.  3] 

3  =  2.02 
and 

■  =0.9 

which  fall  within  the  acceptable  ranges. 

The  initial  plate  used  in  this  study  was  identical  to 
that  used  in  Ref.  3  with  the  exception  that  the  stiffener 
cross  section  used  was  a  tee  rather  than  a  rectangle  (Fig. 

3.2) . 

Since  the  spacing  of  the  stiffeners  in  this  tee- 
stiffened  plate  was  identical  to  that  in  the  original  plate, 


through  application  of  the  kinematic  relationships.  An 
equivalency  can  be  drawn  between  the  motions  of  the  struc¬ 
ture  between  anv  two  instants  in  time  t~  and  t,  and  the 

extrema  of  the  energy  functional  ,'tl(U  -  K)dt  with  resoect 

X0 

to  all  admissible  variations  of  the  strain  energy  U  and  the 
kinetic  energy  K  which  satisfy  the  geometric  boundary  condi¬ 
tions  and  continuity  conditions.  The  equilibrium  equations 
can  then  be  derived  from  the  first  variation  of  the 
energy  functional  through  partial  integration  with  respect 
to  time  and  the  spatial  variables. 

The  variational  approach  provides  a  particularly  effec¬ 
tive  means  for  analyzing  continuous  systems  because  of  the 
way  in  which  it  generates  and  accounts  for  natural  boundary 
conditions  [Ref.  11].  The  disadvant age'  of  the  method  is 
that  it  is  limited  to  conservative  systems.  However,  by 
applying  the  principle  of  virtual  work  STAGS  is  able  to 
greatly  extend  the  classes  of  problems  to  which  it  can  be 
app 1 ied . 

In  the  second  case;  where  inelastic  deformations  such  as 
plastic  strains  or  thermal  expansions  are  encountered,  STAGS 
treats  them  as  "pseudo- loads . "  In  this  manner,  by  applying 
these  pseudo-load  terms  to  the  force  side  of  the  equation, 
an  initial  estimate  of  the  inelastic  deformation  values  can 
be  made  and,  if  convergence  is  obtained,  the  nonconservative 
problem  can  be  solved  by  iteration  upon  a  series  of 


Shock  Analysis  (USA)  is  a  multimodule  routine  which  is 
inserted  as  an  intermediary  and  adjunct  to  the  program  pre¬ 
analyzer  ( STAGS1 )  and  analyzer  ( STAGS2 ) .  The  STAGS1,  STAGS 2 
combination  can  stand  alone  for  some  types  of  general  struc¬ 
tural  analysis.  As  a  result  of  its  add-on  nature,  the  USA 
code  cannot  be  used  in  a  stand  alone  mode,  but  it  is  not 
STAGS  dependent.  Versions  have  been  created  for  use  with 

g 

other  finite  element  structural  analysis  codes. 

B.  THEORY 

The  characteristics  which  make  USA-STAGS  unique  lie  in 
its  applicability  to  a  wide  variety  of  problems  ar.d,  partic¬ 
ularly,  in  the  individual  techniques  used  to  solve  those 
problems.  The  broad,  underlying  theory  behind  those  solu¬ 
tions  is  in  itself  common  to  most  finite  element  codes. 

Briefly,  the  technique  used  by  codes  of  this  genre  is  to 
solve  the  differential  equations  of  motion  for  a  structure 
in  order  to  determine  its  responses  to  its  environment 
[Ref.  10:  pp.  2-1  through  2-14].  By  then  applying  the 
constitutive  equations  (Hooke's  law  generalized  to  include 
anisotropic  materials),  the  stresses  in  the  governing  equa¬ 
tions  of  motion  may  be  expressed  in  terms  of  strains.  These 
strains  may  in  turn  be  expressed  as  the  first  partials  of 
the  displacements  with  respect  to  the  spatial  coordinates 

g 

As  of  this  writing,  USA  has  been  linked  to  N’ASTRAN , 
GENSAM,  SPAR,  and  DIAL. 
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IV.  THE  USA-STAGS  CODE 


A .  GENERAL 

Structural  Analysis  of  General  Shells  (STAGS)  is  the 
finite  element  code  first  developed  by  Lockheed  Missiles  and 
Space  Company,  Inc.  during  the  1960's  in  response  to  an 
initial  tasking  by  the  Navy  and  the  Air  Force  plus  an  inter¬ 
nally  perceived  need  for  a  general  structural  analysis 
program . 

The  current  version  of  STAGS  is  STAGS  C-l ,  a  dynamic 
analysis  code  with  broad  applications  to  the  aerospace 
industry  which  has  also  found  application  in  ship  design  and 
engineering.  A  new  version,  STAGS  2000,  is  now  under  develop 
ment  which  will  have  an  entirely  changed  data  manager 
concept  from  that  used  in  the  C-l  version.  This  change 
should  make  the  new  code  both  easier  to  use  and  far  more 
economical  of  the  user's  time.  STAGS  2000  will  appear 
initially  as  a  static  analysis  code,  but  will  eventually 
exceed  STAGS  C-l's  capabilities. 

STAGS  C-l  is  supported  by  Computer  Software  Management 
and  Information  Center  (COSMIC)  under  the  sponsorship  of  the 
National  Aeronautics  and  Space  Administration  (NASA). 


The  USA-STAGS  code  is  an  extension  of  the  basic  STAGS 
routine  to  permit  evaluation  of  impulsive  loadings  upon 
design  structures  in  a  fluid  medium.  As  such,  Underwater 


closer  investigation  of  the  shock  effects  themselves, 
without  attempting  either  to  retain  a  set  plastic  modulus 
or  to  approximate  ships'  grillage. 

The  result  of  this  simplification  was  a  thickening  of 

the  base  panel  to  i  inch  to  increase  resistance  to  shear. 

The  stiffener  was  restored  to  the  original  rectangular 

configuration  but  the  thickness  was  increased  to  J  inch 

7 

and  the  depth  to  1J  inches.  Most  important,  in  this 
attempt,  a  single,  longitudinally  oriented  stiffener  was 
used  rather  than  the  two  transverse  stiffeners  in  an  effort 
to  provide  greater  flexibility  along  the  length  of  the 
stiffener  (Fig.  3.6). 

Although  the  charge  used  to  test  the  rectangular  stif¬ 
fened  plate  was  small  (£  pound  of  TNT),  and  the  correspond¬ 
ing  plate  deflections  were  far  too  slight  for  any  firm 
conclusions  to  be  drawn,  it  can  be  seen  in  Figures  3.7  and 
3.8  that  tripping  action  appears  to  have  begun  to  manifest 
itself.  The  single,  longitudinally  oriented  stiffener  seems 
to  be  the  superior  model;  and,  while  further  testing  is 
required  at  larger  charge  weights,  its  employment  is  tenta¬ 
tively  recommended  for  future  studies. 

7 

This  decision  was  made  somewhat  arbitrarily.  Dimen¬ 
sions  were  chosen,  in  part,  to  reduce  the  cost  and  simplify 
machining  the  plate.  It  is  interesting  to  note,  however, 
that,  with  3  =  1.52  and  \  =  1.1,  this  plate  is  not  far 
outside  the  acceptable  grillage  range. 


is  applied  to  the  top  fibers  of  the.  stiffener.  Recognizing 

that  a  flat  plate  will  tend  to  deform  away  from  the  source 

of  an  explosion,  the  plates  in  this  series  of  studies  have 

been  oriented  with  the  stiffener  sides  outward  into  the 
5 

fluid  medium. 

As  stated  earlier  in  this  chapter,  the  original  intent 

in  this  series  of  studies  was  to  maintain  an  essentially 

fixed  plate  model,  varying  only  the  stiffener  geometries. 

Unfortunately,  the  results  of  the  first  two  explosion  tests 

(on  the  original  rectangular-st if f ened  plate  [Ref.  3],  and 

on  the  tee-stiffened  plate  here)  were  not  satisfactory  in 

view  of  the  stated  objectives.  The  two  transversely  mounted 

stiffeners  were  too  short  and  the  entire  assembly  too  rigid 

to  effectively  demonstrate  the  desired  stiffener  tripping. 

The  shock  wave  essentially  "punched"  through  the  panel 

material,  shearing  it  at  the  fixed  boundaries  and  at  the 

0 

base  of  the  stiffeners  (Fig.  3.5). 

The  results  of  the  two  initial  tests  indicated  that  a 
design  change  was  in  order.  The  third  plate  of  the  series 
(the  second  of  this  study)  was  a  simplified  model  to  permit 

5 

This  is  hardly  the  expected  orientation  on  a  surface 
vessel,  but  is  common  on  many  submarine  pressure  hulls. 

6 

Although  not  easy  to  see  in  Figure  3.5,  not  all  of  the 
damage  to  the  tee-stiffened  plate  was  done  by  the  incident 
shock  wave.  Apparently,  when  the  bottom  reflection  hit  the 
plate  and  its  backing  structure,  it  was  forced  upward  into 
its  own  supporting  cables.  Much  of  the  damage  to  the  stif- 
finers  in  particular  seems  to  have  been  caused  by  the 
actions  of  that  cable. 


which  gives  a  radius  of  curvature: 


=  [0.092562  in^ / 1.0256  in2]* 

=  0.300419  in 

Thus,  for  the  tee-stiffened  plate, 

=  12  40000  * 

"  -300419TT  1Q  x  1q6 

=  0.804144 

which  is  well  within  the  parameter  range  required  for 
modelling  ships'  grillage. 

The  second  requirement  levied  upon  these  test  plates  was 
the  ability  to  demonstrate  stiffener  tripping. 

Stiffener  tripping  has  been  identified  as  the  lateral- 

« 

torsional  instability  of  the  stiffener  as  it  becomes 
suddenly  unstable  and  fails  under  load.  It  is  identified  by 
a  characteristic  warping  and  buckling  of  the  stiffener.  It 
is  significant  that  the  continued  resistance  of  a  stiffened 
structure  to  deformation  under  loading  is  dramatically 
reduced  after  the  stiffeners  have  tripped.  In  earlier 
studies  of  impulsive  loadings  upon  ships'  grillage,  tripping 
has  been  identified  as  "...  a  primary  ductile  failure  mode 
for  ship  structure."  [Ref.  9:  p.  2]. 

Previous  observations  have  indicated  that  tripping  is 
most  readily  initiated  when  an  axially  compressive  loading 
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=  0.607747  in. 


Plascic  Section  Modulus: 


Z  =  Sum  of  the  moments  about  the 
P  neutral  axis  of  the  cross 

sections  above  and  below  the 
axis 


-)l  K  ^  i ,  'o  j 


Zp  =  (bjh^  (h2  +  ^h  1 ) 

h^  *  1/8  in-  =  0.125000  in. 
bx  -  in.  =  0.718750  in. 

h2  -  ||  -  0.607747 

=  0.173503  in. 


h3  »  0.607747  in. 

b„  =  b,  =  ^  in.  =  0.125000  in. 
2  Jo 


(0.718750) (0. 125) (0. 173503  +  4  -  0.125) 
+  4  (0. 125) (0. 173503)2 
+  4  (0 . 125) (0 . 607747 ) 2 
0.046170  in3 


clearly  the  3  values  remained  invariant.  Given  the 
changed  geometry  of  the  stiffener  the  effective  radius  of 
gyration  of  the  plate /st if fener  assembly  (k)  could  be 
expected  to  change,  and  with  it  the  value  of  X. 

Before  that  value  X  could  be  determined,  however,  the 
dimensions  of  the  stiffener  itself  first  had  to  be 
established . 

In  an  effort  to  ensure  that  stiffener  responses  were  a 

function  of  stiffener  design  alone,  the  plastic  moduli,  Z  , 

of  the  two  stiffener  cross  sections  were  set  as  close' as 

possible  to  equal.  For  the  original,  rectangular-st if fened 

3 

plate,  Zp  equaled  0.046875  in  . 

By  a  process  of  iteration,  a  tee-stiffener  with  flange 
dimensions  23/32  inch  by  1/8  inch,  and  web  dimensions  25/32 
inch  by  1/8  inch  was  fixed  upon  as  giving  an  elastic  modulus 

3 

Z  =  0.046170  in  ,  which  is  within  1.5%  of  that  for  the 

p 

rectangular  stiffener  (Fig.  3.3). 

Now,  applying  these  dimensions  to  the  grillage  param¬ 
eter,  X,  we  find  that  the  neutral  axis  of  the  effective 
stiffener  (Fig.  3.4)  is: 


=  0.2220  in 


The  moment  of  inertia  of  the  effect ive  stiffener  is: 

2 
4 


I  =  ZI.  +  ZA.d. 
ZZ  i  11 


=  0.092562  in 


conservative  problems  [Ref.  10:  pp .  2-7,  2-8].  The  range 
of  cases  to  which  STAGS  may  now  be  applied  is  shown  in 
Figure  4.1. 

After  being  reduced  to  functions  of  displacement  only, 
the  governing  equations  of  motion  for  a  discretized  system 
may  be  expressed  in  the  form: 

MsH  +  Ds(X)  +  BS(X>  +  Kg(x)  =  f 

where  x  is  the  vector  of  displacement  components,  M  is 
~  s 

the  structural  mass  matrix,  f  is  the  vector  of  external 
forces,  and  Kg  is  the  (generally  nonlinear)  structural 
stiffness  operator.  The  operators  Bs  and  Dg  include 
forces  that  are  functions  of  structural  deformation  and 
deformation  velocity,  respectively. 

Equation  4.1  is  the  matrix  equation  that  STAGS  must 
solve  to  find  structural  displacements.  For  fluid/struc¬ 
tural  interactions  the  generalized  force  vector  becomes  a 
much  more  complex  entity  [Ref.  12:  pp .  2-1  through  2-12]. 

For  submerged  structures  excited  by  an  acoustic  wave, 
the  exciting  forces  are  given  by: 

f  =  -GA^(Pj  +  Pg)  (eqn  4.2) 

where  Pj  and  Pg  are  the  incident  and  scattered  nodal  pres¬ 
sure  vectors  acting  upon  the  surface  of  the  structure  at  the 
fluid-structural  interface.  A ^  is  the  diagonal  area  matrix 
associated  with  elements  in  the  fluid  mesh,  and  G  is  the 
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matrix  which  relates  the  structural  and  fluid  surface 
forces.  It  can  be  seen  from  the  development  of  the  pressure 
characteristics  of  an  underwater  shock  wave  in  Chapter  1, 
that  for  an  underwater  explosion  the  incident  pressure 
vector  Pj  is  known.  The  scattered  pressure  vector  Pg  on  the 
other  hand  is  dependent  upon  the  configuration  of  the  struc¬ 
tural  surface  and  is  therefore  not  only  not  immediately 
known  but  is  also  changing  as  the  structure  deforms. 

It  is  to  solve  this  computational  dilemma  and  thereby 
solve  the  combined  equation  of  motion: 

M  x  +  Dg(x)  +  Bg(x)  +  Kg(x)  -  -GAf(Pg  +  ?S}  (eqn  4*3) 

that  the  USA  code  was  created. 

The  USA  code  finds  Pg  by  utilizing  the  Doubly  Asymptotic 
Approximation  ( DAA )  and  a  staggered  solution  procedure. 

The  DAA  is  an  approximate  relation  which  approaches 
exactness  at  both  the  high-frequency  (early-time)  and  low- 
frequency  (late-time)  limits  to  the  passage  of  the  shock 
front.  It  may  be  written: 

MfPg  +  pcAfPg  =  pcMfVg  (eqn  4.4) 

where  vg  is  the  vector  of  scattered-wave  fluid-particle 
velocities  normal  to  the  structural  surface  at  the  fluid- 
structural  interface.  is  the  fluid  mass  matrix  for  the 

fluid  mesh  at  the  interface,  and  p  and  c  are  the  density  and 
sound  velocity  of  the  fluid,  respectively. 


When  fluid  and  structure  are  excited  by  an  acoustic 
wave,  the  incident  and  scattered  fluid  particle  velocities 
(Vj  and  Vg)  may  be  related  to  the  structural  response  by  the 
kinematic  compatability  relation: 

T  • 

G  x  =  Vj  +  Vg  (eqn  4.5) 

Solving  for  vg  and  substituting  the  resultant  into  equation 
4.4  yields 

MfPg  +  pcAfPg  =  ccMf(GTx  -  Vj)  (eqn  4.6) 

Equations  4.3  and  4.6  are  now  each  expressed  in  terms  of  the 
single  unknown  Pg  and  a  solution  may  be  obtained  by  a  stag¬ 
gered  iteration  process  between  the  two  equations  for  each 
time  step. 

C.  ORGANIZATION 

As  has  already  been  stated,  USA-STAGS  is  a  modular  code. 
Functionally,  it  may  be  divided  into  five  basic  components. 

1.  The  structural  preprocessor 

2.  The  fluid  mass  preprocessor 

3.  The  augmented  matrix  preprocessor 

4.  The  time  integration  processor 

5.  The  data  postprocessor 

Each  of  these  components  is  run  individually  in  sequence 
utilizing  the  outputs  from  the  previous  modules. 


1 .  The  Structural  Preprocessor 

The  structural  preprocessor  is  the  unit  that 
assembles  the  structural  mass  and  stiffness  matrices, 
creates  the  finite  element  mesh,  and  otherwise  provides  an 
encoded  description  of  the  structure's  internal  and  ex¬ 
ternal  degrees  of  freedom  [Ref.  12:  pp .  3-3,  3-4]. 

For  USA-STAGS,  the  structural  preprocessor  is  the 
STAGS1  portion  of  the  STAGS  code  itself. 

Input  for  STAGS  is  made  to  each  of  ten  functional 
subdivisions  which  correspond  to  the  types  of  input  that  go 
into  each.  These  subdivisions  are: 

-  Summary  and  Control  Parameters 

-  Computational  Strategy  Parameters 

-  Data  Tables 

-  Geometry  and  Material 

-  Discretization 

-  Discrete  Stiffeners* 

-  Boundary  Conditions 

-  Loads* 

-  Output  Control 

-  Finite  Element  Units* 

(Subdirectories  designated  with  *  are  not  required  with  all 
models . ) 

Input  to  each  subdirectory  is  in  the  form  of  a 
series  of  data  records  consisting  of  real  and  integer 
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numerical  data  :  : - 


discussed  in  greater 


detail  in  t  he  f 

2  .  The  FI  1 :  i  T -  :  •  ■  :_L__  r 

The  f  1  uu  ;  T.a.-  •  r  •  r  rI.LT.IA5  '  constructs  the 

fluid  mass  matrix  f.r  a  ^*ru:*ure  submerged  in  a  fluid. 
Inputs  include  [Ref.  12:  pp.  3-4,  3-5]  the  global  coordi¬ 
nates  for  the  structural  geometry  >  input  from  the  structural 
preprocessor),  plus 

-  Fluid  wet-surface  mesh 

-  Element  definitions 

-  Material  properties  of  the  fluid 

-  Location  of  free  surface 

-  Model  symmetries 

3  .  The  Augmented  Matrix  Prepr-V'ess-  _r 

The  augmented  matrix  prepr-'cess-.-r  A'.'GMAT  receives 
the  output  from  the  structural  and  f’.ui  i  ma->>  r  'e^^rs  t  :> 
create  the  matrices  required  for  so  1  u*  i  -r.  :  ■  • t  jgmented 

equations  4.3  and  4.6. 

4.  The  Time  Integration  Pro r 

The  funct  ional  heart  of  t  he  ' '  5  A  • .  4  .• 

time  integrator  processor  (TIMIN'!  1  ms.-.  :  .  *  •  ■  m.  :  pig¬ 
ments  the  structural  processor  STAGS2  to  imp  l-men *  *  he  stag¬ 

gered  solution  technique;  creating  the  solution  data  files 
over  the  time  range  specified.  Input  includes  the  data 
passed  from  AUGMAT  plus 


-  Spatial  geometry  of  charge  and  target 

-  Initial  incident  pressure  and  pressure  wave  decay 
characteristics 

-  Time  integration  information 

-  Restart  data 

TIMINT  includes  an  internal  postprocessing  capability. 
Generally,  however,  this  function  is  left  to  the  following 
module . 


5.  The  Data  Postprocessor 


USA-STAGS  Provides  a  selection  of  three  possible 


postprocessors:  POSTP ,  POSTPR ,  and  STAPL .  Their  uses  will 


be  discussed  in  the  following  section. 


D.  USING  THE  CODE 

As  has  already  been  stated,  the  STAGS  code  is  extensive 
in  its  concept;  attempting  tp  cover  in  a  single  code  as 
many  approaches  to  problem  solution  as  can  be  linked  in 
a  single  contextual  framework. 

It  is  beyond  the  scope  of  this  work  to  cover  all 
possible  aspects  of  the  code;  indeed  within  the  stated  capa¬ 
bilities  of  STAGS  this  approach  may  seem  peripheral,  being 
limited  as  it  is  to  those  aspects  clearly  related  to  under¬ 
water  shock  analysis  of  flat  plate  models. 

STAGS  is  designed  to  be  used  for  any  one  of  six  types  of 
analysis  [Ref.  13:  p.  2-1]. 

-  Linear  stress  analysis 

-  Geometrically  nonlinear  elastic  stress  analysis 


-  Inelastic  stress  analysis,  geometrically  Linear  or  • 
nonlinear 

-  Bifurcation  buckling  analysis  with  linear  or  geometri¬ 
cally  nonlinear  stress  state  (elastic) 

-  Small  vibration  analysis  based  on  linear  or  geometri¬ 
cally  nonlinear  stress  state  (elastic) 

-  Transient  response  analysis,  linear  or  geometrically 
nonlinear,  elastic  or  inelastic 

In  this  study,  only  two  of  those  capabilities  were  used: 

(a)  Small  vibrations  (in  a  stress  free  state),  and  (b) 
Geometrically  nonlinear  transient  response  analysis.  Of 
these,  the  second  alone  is  truly  applicable  to  USA-STAGS . 

The  small  vibration  analysis  was  used  as  a  check  on  the  code 
and  input  data.  Since  the  actual  inputs  for  the  two  anal¬ 
yses  vary  but  slightly,  both  shall  be  addressed  together. 

Before  that  can  be  done,  however,  the  codes  must  somehow 
be  brought  into  some  semblance  of  conformity;  for  they  are 
not  normally  provided  as  a  single  unit. 

Both  codes  are  provided  in  tape  form:  STAGS  from 
COSMIC,  and  USA  directly  from  Lockheed.  These  are  designed 
to  be  loaded  onto  the  VAX  11/780  computer  using  the  standard 
MOUNT  MF :  command . 9 

It  is  suggested  that  ,  at  least  during  the  early  stages 
of  processing,  each  tape  be  assigned  to  a  separate  direc¬ 
tory,  for  many  of  the  source  files  will  be  culled  as  the 


executable  files  (*.EXE)1(1  are  created.  The  source  files 
supplied  with  each  code  are  specific  to  the  code  until  the 
final  linkage,  and  their  sneer  number  can  cause  confusion  if 
not  approached  methodically. 

The  STAGS  tape  as  provided  contains  all  the  necessary 
ingredients  for  the  functioning  codes  except  the  executables 
themselves.  These  are  left  to  be  linked  on  the  individual 
computer.  The  link  commands  for  creation  of  these  executa¬ 
bles  are  provided  with  the  tape  on  the  INFO. VAX  file,  which 
also  contains  a  full  listing  of  the  tape's  contents.  The 
commands  that  are  of  immediate  interest  are  LSTAGS1.COM, 
LSTAGS2.COM,  LPOSTP.COM,  and  LSTAPL.COM.  These  should  be 
copied  into  the  file  and  executed  in  the  normal  manner 
for  a  program  run  from  the  terminal.11  The  result  will  be 
the  executable  files  for  the  respective  module. 

Executables  are  provided  with  the  USA  tape;  however, 

since  relinking  may  be  required  to  adapt  the  code  to  local 

12 

filenames,  link  command  files  are  also  provided.  The 

1(10n  the  VAX,  the  asterisk  represents  a  "wildcard”  which 
may  be  substituted  for  any  filename,  filetype,  or  version 
number  (the  complete  listing  for  a  file  is:  filename, 
filetype,  version  number)  whenever  a  general  category  is 
desired  or  the  listing  is  to  be  repeated. 

11This  is  done  by  simply  typing  ”0  command  filename"  and 
hitting  the  return  key. 

12 

The  USA  tape  also  includes  another  particularly  nice 
feature,  the  complete  user's  manual  for  each  of  the  modules. 
This  ensures  that  the  most  current  version  is  available; 
and,  more  important,  the  version  of  the  manual  that  is 
pertinent  to  that  tape.  A  vital  accompaniment  since  the 
tapes  are  continually  being  updated. 


command  files  to  be  used  in  constructing  the  USA  executables 


are  LINKFLU.COM,  LINKAUG.COM,  LINKUSAS.COM,  and  LINKPOS.COM. 
It  is  the  LINKUSAS.COM  file  which  actually  performs  the 
marriage  between  the  codes;  allowing  the  integration/ 
structural  analysis  functions  of  STAGS  to  be  coupled  with 
the  USA  fluid  processor. 

Once  all  of  the  executable  files  have  been  created  and 
tested  in  a  full  run  of  the  combined  code;  all  files  but 
the  executables  may  be  deleted,  thus  realizing  a  considerable 
increase  in  computational  space.  ^ 

The  basic  strategy  that  will  be  followed  here  shall  be 
to  address  USA-STAGS  as  a  single  code  broken  into  five 
different  modules  (see  preceding  section);  each  of  which  can 
be  run  independently  in  stepping  stone  fashion  once  neces¬ 
sary  input  from  earlier  modules  has  been  obtained. 

The  first  module  to  be  run  is  STAGS1 .  The  actual  input 
for  this  includes  all  information  needed  for  both  STAGS1  and 
STAGS2  execution  in  a  single  entry  file. 

The  first  five  data  records:  A-l,  B-l,  B-2  ,  B-3,  and 
C-l  are  common  to  all  STAGS  runs.  The  inputs  as  set  forth 
in  the  COSMIC  STAGS  users  guide  [Ref.  13:  pp.  3-1  through 
3-11]  are  largely  self-explanatory  but  a  few  points  should 
be  noted. 


13 

Total  space  required  for  the  executables 
than  4000  blocks;  compared  to  the  almost  13000 
the  executables  plus  source  codes. 


alone  is  less 
blocks  for 
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First,  IPOSTO ,  IPOST1 ,  and  IPOST2  are  applicable  only  to 
the  STAGS  postprocessors  POSTP ,  and  STAPL .  If  the  user  is 
interested  only  in  the  nodal  displacements  resulting  from 
underwater  shock  these  may  be  derived  directly  from  either 
of  the  USA  modules  TIMINT  or  POSTPR.  Thus  IPOSTO,  IPOST1, 
and  IP0ST2  may  be  set  to  zero  thereby  suppressing  creation 
of  model  and  solution  data  files  ( * . MOD  and  F0R022.DAT)  and 
realizing  some  savings  in  computer  space.  It  is  suggested 
that  this  be  done  in  any  case  during  early  model  development 
and  while  testing  code  inputs  (as  with  the  small  vibrations 
analysis).  If,  however,  strains  are  required;  these  are 
computed  only  within  the  STAGS  postprocessors  and  IPOSTO 
must  be  set  to  1.  Considerable  amounts  of  computer  time  and 
space  can  be  used  in  processing  a  large  model;  time  which 
would  be  wasted  if  this  is  not  done.  Similarly,  IP0ST1  and 
IP0ST2  should  be  set  to  a  non-zero  value  if  displacement  and 

stress  plots  are  desired  at  the  IPOSTlst  and  IPOST2nd  load 

*  14 

steps . 

Second,  it  should  also  be  noted  well  that  IFLU  must  be 
set  to  1  if  underwater  shock  analysis  using  USA  is  to  be 
undertaken.  This  will  cause  the  FOR003.DAT  and  FOR004.DAT 
files  required  for  USA  to  be  created.  These  files  are  not 


14 

This  is  required  only  when  using  STAPL.  POSTP  is 
capable  of  calculating  displacements  and  stresses 
internally . 
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unduly  lengthy,  but  for  the  sake  of  economy  should  be 

suppressed  for  non-USA  runs. 

The  C-l  data  record  is  meaningless  for  small  vibration 

analysis  and  transient  response  analysis  but  must  be 

i  5 

included.  A  value  of  0.  should  be  entered. 

The  only  true  variation  in  data  record  sequence  between 
small  vibration  analysis  and  transient  response  analysis 
comes  in  the  next  group  of  cards.  The  proper  sequencing  for 
small  vibrations  is  D-2 ,  D-3,  F-l;  while  for  transient  anal¬ 
ysis  it  is  D-l,  E-l,  E-2 >  F-l.  These  are  the  Strategy 
Parameters  and  the  difference  lies  in  the  need  for  integra¬ 
tion  over  time  in  the  transient  model.  The  D-l  record  sets 
tolerances,  limits  the  amount  of  CPU  time,  allows  restarts, 

etc.;  while  the  E  records  control  the  actual  time  integra- 
16 

tion  itself. 
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This  must  be  a  real  number.  STAGS  is  often  blind  to 
data  not  submitted  in  the  proper  format.  It  will  not  record 
an  error,  but  will  simply  pass  over  the  improperly  input 
data  until  finding  some  in  the  format  it  seeks.  This  can 
lead  to  sometimes  dangerous  distortions  of  the  input 
records . 

16 

At  the  time  that  this  study  was  initiated  it  was 
understood  that  implicit  integration  using  K.  C.  Park's 
formula  ( IMPL  =  0,  METHOD  =  4  on  the  E-2  record)  was  the 
only  form  of  integration  that  should  be  used  with  USA.  As 
of  this  writing,  however,  Lockheed  informs  that  USA  may  now 
also  be  used  with  explicit  integration.  It  has  been 
suggested  that  this  might  provide  a  faster  and  more  econom¬ 
ical  solution.  Due  to  time  limitations  that  has  not  been 
tested  here. 
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The  D-2  and  D-3  records  are  specific  to  modal  analysis 
and  deal  with  the  number  and  range  of  eigenvalues  to  be 
calculated . 

Data  records  F  through  P  pertain  to  the  model  itself. 

The  records  F  through  L  provide  information  general  to  the 
entire  model  and  are  listed  in  normal  alphabetical  sequence 
once  through  only,  although  specific  cards  may  be  repeated 
within  the  sequence.  M  through  P,  plus  Q  and  R  are  specific 
to  the  individual  shell  and  beam  components  and  are  repeated 
in  sequence  once  for  each  component  (see  examples  in  the 
appendix ) . 

The  F-l  and  F-2  records  provide  grid  point  and  stiffener 
summaries,  respectively.  The  information  in  F-l  controls 
generation  of  the  finite  element  mesh.  Inputs  are  simply 
the  number  of  nodal  rows  and  columns  for  each  shell.  Actual 
choice  of  which  shell  direction  shall  consist  of  rows  and 
which  columns  is  unimportant  as  long  as  that  choice  remains 
consistent  for  the  local  X  and  Y  coordinates  across  the 
entire  model  (rows  are  lines  of  consistent  X,  columns  are 
lines  of  consistent  Y).  Entry  is  made  sequentially;  number 
of  rows(I),  number  of  columns(I)  for  each  shell  in  order  of 
input.  The  mesh  created  forms  the  master  elements  which 
will  then  be  fitted  into  the  structural  geometry  by  the 
information  contained  in  the  M  records. 

The  F-2  record  is  a  summary  of  shell  stiffeners.  Use  of 
this  record  could  be  extremely  useful  in  determining 
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generalized  responses  on  large  scale  models,  but  was  incon¬ 
sistent  with  the  purposes  of  this  study  because  of  the 
inherent  assumption  that  stiffeners  do  not  deform. 

Record  G-l  establishes  the  compatibility  conditions  for 
the  various  shells.  The  instructions  listed  in  the 
COSMIC-STAGS  user's  manual  are  straightforward,  but  refer¬ 
ence  to  Appendix  B  may  prove  helpful  for  establishing  proper 
linkages.  It  will  be  noted  that  in  each  case,  for  a  flat 
plate,  side  1  is  that  nearest  and  parallel  to  the  global  Y 
axis  while  side  4  has  a  similar  correspondence  to  the  global 
X  axis. 

Records  I,  J,  K,  and  L  constitute  the  Data  Table  section 
of  STAGS.  Again  the  user’s  manual  is  self-explanatory  and 
will  not  be  explicated  except  to  note  that  the  Material 
Record,  I — 1 ;  the  Material  Property  Record,  1-2;  and  the 
Plasticity  Record  1-3  (if  used);  are  listed  in  sequence 
once  for  each  of  the  materials  enumerated  in  the  B-3  record. 

If  plastic  behavior  is  expected  in  the  model  and  the 
White-Bessel ing  theory  of  plasticity  is  chosen  ( IPLST 
options  0  or  1  on  the  1-1  record)  especial  note  should  be 
made  of  the  user's  manual's  warnings  concerning  the  slope  of 
the  line  connecting  points  submitted  for  the  discretized 
stress-strain  curve  in  the  1-3  record.  The  change  in  slope 
at  each  point  must  be  negative.  Another  point  not  mentioned 
in  the  user's  manual,  but  which  should  definitely  be  noted, 


is  that  the  first  stress-strain  data  entry  on  the  1-3  record 


necessitated  by  the  use  of  a  half-plate  model  instead  of  a 
quarter-plate  set  the  total  number  of  fluid  elements  beyond 
i hat  number  when  the  0.6  inch  element  was  used.  Rectangular 
elements  1.2000  inches  by  1.2857  inches  were  chosen  for 
the  base  shells  in  order  to  allow  strain  gage  placement  at 
the  center  of  each  shell.  Elements  for  the  stiffener  shell 
*ere  1.2857  inches  by  0.6250  inch.  This  corresponded  to  6 
nodal  rows  and  8  nodal  columns  for  shell  units  1  and  2, 
and  3  rows  and  8  columns  for  shell  unit  3  (Fig.  4.5). 


Figure  4.5  Half  Model  of  Rectangular-St if fened  Plate 
For  the  final  model  of  both  the  tee-stiffened  and 
rectangular- stiffened  plates,  the  finite  element  mesh  was 
placed  at  the  mid-planes  of  the  shell  units.  This  was  not 


81 


extensive  rewriting  could  be  avoided  throughout  the 

22 

USA-STAGS  modules. 

The  STAGS  code  will  generate  the  finite  element  mesh. 

The  user  need  only  specify  the  number  of  nodal  rows  and 

columns.  In  the  case  of  the  tee  stiffened  plate,  an  element 

0.6  inches  by  0.6  inches  was  chosen  in  order  to  keep  the 

mesh  as  fine  as  possible  and  yet  be  able  to  place  the  active 

filaments  of  the  strain  gages  used  during  the  underwater 

23 

explosion  testing  entirely  within  the  element.  This  corre 
sponded  to  11  nodal  rows  and  6  nodal  columns  in  shell 
unit  1,  11  rows  and  11  columns  in  shell  unit  2,  and  11 
rows  and  3  columns  in  each  of  the  three  stiffener  shells 
(Figure  4.3). 

Since  the  stiffener  used  on  the  second  test  plate  was 
rectangular,  three  shell  units  sufficed  to  model  it  (Fig. 
4.4).  At  the  time  that  this  plate  stiffener  was  discre¬ 
tized,  the  limit  set  upon  the  number  fluid  mesh  elements  by 

24 

the  FLUMAS  module  was  171.  The  increased  surface  area 

22 

Unfortunately,  it  was  not  possible  to  carry  through 
with  this  plan  since  the  plate  with  the  longitudinal  rectan¬ 
gular  stiffener  completely  changed  material  dimensions  and 
permitted  only  one  axis  of  symmetry. 

23 

STAGS  computes  stresses  and  strains  at  the  geometric 
centroids  of  the  elements.  Strain  gages  should  therefore  be 
centered  on  the  element . 

24 

It  turns  out  that  this  is  relatively  easy  to  change, 
but  it  involves  technical  considerations  that  are  beyond  the 
scope  of  this  work.  Questions  about  such  modification 
should  be  directed  to  Lockheed.  The  limit  on  most  accounts 
has  recently  been  increased  to  400. 


E.  DISCRETIZING  THE  MODEL 


As  its  name  implies,  the  primary  building  unit  in  the 
STAGS  code  is  the  shell  element.  These  elements  may  be 
constructed  into  shapes  which  are  warped  through  almost  any 
conceivable  configuration  in  three-space  and  beam  elements 
may  be  added,  but  the  basic  unit  remains  two-dimensional. 

The  flat  plate  models  constructed  for  this  study,  then, 
were  ideally  suited  fo^  discretisation  using  such  elements. 

Each  plate  with  its  associated  stiffeners  was  capable  of 
being  divided  into  several  flat  shell  units  composed  of 
square  or  nearly  square  elements.  Additionally,  the  even 
spacing  of  the  stiffeners  allowed  advantage  to  be  taken  of 
the  plate  symmetries  to  reduce  the  overall  size  of  the  model 
and  corresponding  length  of  the  computer  run. 

The  transverse  tee-stiffened  plate  established  symmetry 
in  two  directions  permitting  a  quarter  model  to  be  used. 

That  quarter  model  was  then  divided  into  five  shells;  two 
for  the  base  plate,  and  three  for  the  stiffener.  These 
shells  were  numbered  beginning  with  the  two  shells 
comprising  the  base  upward  through  the  stiffener  (Fig. 

4.2).  This  was  done  to  permit  changes  to  be  made  to  the 
code  input  by  simply  adding  or  deleting  shells.  As  long  as 
the  basic  plate  dimensions  and  symmetries  remained  the  same, 


and  may  be  renamed  to  the  STAGS  subdirectory  for  use  in 
either  of  these. ^ 

Both  POSTP  and  STAPL  have  been  designed  to  provide  not 

only  displacements  but  also  stresses  and  strains;  POSTP  in 

tabular  form  and  STAPL  in  graphical  on  a  CALCOMP  plotter. 

At  the  time  of  this  writing,  the  Naval  Postgraduate  School 

did  not  have  a  graphical  capability  on  VAX/VMS  so  neither 

the  POSTPR-graphical  nor  STAPL  options  were  run  at  the 
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Naval  Postgraduate  School. 

In  accordance  with  common  usage  on  the  VAX,  the  sequence 
of  commands  required  to  execute  any  module  have  been  entered 
into  command  files  (*.COM).  Examples  are  given  in  Appendix 
A. 

Once  the  appropriate  file  names  have  been'  entered  into 
the  command  file  for  the  module  to  be  run ,  the  command  file 
may  be  executed  by  entering  command  filename"  for  runs 
from  the  terminal;  or  "SUBMIT  filename"  for  batch  runs.  A 
file  definition  statement  (SET  DEF  [*.*])  is  required  for 
batch  runs  but  is  not  necessary  for  runs  from  the  terminal. 

20 

It  should  be  noted  that  the  TIMINT  entry,  NSAVER,  is 
the  frequency  with  which  responses  are  saved  for  the  USA 
postprocessor  only.  The  frequency  with  which  data  is  saved 
for  the  POSTP  and  STAPL  postprocessors  is  actually  tied  to 
the  restart  data  in  NRESET.  If  stress  and  strain  informa¬ 
tion  is  required  at  each  timestep,  NRESET  should  therefore 
be  set  to  1. 

2 1 

POSTPR  and  STAPL  plots  shown  in  later  sections  of  this 
work  were  done  by  Lockheed. 


FLUMAS  input  must  be  changed  accordingly  if  a  vertical 
aspect  is  required. 

Provision  has  also  been  made  in  TIMINT  to  print  nodal 
displacements  without  recourse  to  a  post-processor.  It  is 
recommended  that  this  option  not  be  exercised.  A  run  for 
any  of  the  other  USA-STAGS  modules  can  be  made  in  a  matter 
of  minutes.  Running  TIMINT  can  easily  consume  several  hours 
of  CPU  time  for  a  moderately  sized  model.  It  is  therefore 
not  time  or  cost  effective  to  run  TIMINT  each  time  a  new  set 
of  displacements  is  needed. 

The  solution  data  file,  POSNAM,  output  by  TIMINT  is 

specific  to  the  USA  post-processor,  POSTPR.  If  nodal 

displacements  only  are  desired  this  is  by  far  the  easiest 

and  most  economical  route  to  follow.  It  is  not,  however, 

the  only  option.  USA-STAGS  provides  a  selection  of  three 

possible  post-processors.  As  mentioned,  POSTPR  is  internal 

to  USA  and  may  be  used  to  provide  either  data  listings  or 

graphical  displays  of  nodal  displacements.  The  other  two, 

POST?  and  STAPL,  are  a  part  of  the  STAGS  package.  A  second 

19 

solution  data  file  in  STAGS  format,  FOR022.DAT,  is  created 

19 

It  is  actually  created  by  STAGS2  during  the  TIMINT 

run . 


The  module  AUG.MAT  is  the  component  of  USA  which  assem¬ 
bles  the  data  from  the  FLUMAS  module  and  from  STAGS1  into 
the  specific  constants  and  arrays  required  for  the  staggered 
solution  process.  The  input  deck  for  this  module  is  very 
short  and  experience  indicates  that  problems  with  the  module 
will  generally  be  caused  by  input  errors  to  either  STAGS1  or 
FLUMAS . 

USA's  time  integrator  module,  TIMINT,  links  the  results 
of  the  AUGMAT  module  which  are  stored  in  the  data  file 
PRENAM  back  to  STAGS2  and  then  unites  with  STAGS2  to  conduct 
a  step-by-step  numerical  time  integration  of  the  governing 
equations  for  submerged  structures  exposed  to  shock  waves. 
The  time  integration  information  input  to  TIMINT  will  over¬ 
write  earlier  entries  to.  the  STAGS  deck.  It  is  suggested, 
however,  that  these  entries  should  be  made  consistent. 

Input  to  TIMINT  includes  not  only  time  integration 
information  but  also  the  spatial  coordinates  for  the  charge 
and  for  the  point  on  the  structure  closest  to  the  charge 
(the  standoff).  These  distances  are  made  in  reference  to 
the  origin  of  the  global  structural  coordinate  as  defined  in 
the  STAGS  input.  One  point  of  caution:  The  surface-cutoff 
provision  in  FLUMAS  assumes  a  horizontal  separation  between 
charge  and  target.  As  the  charge  approaches  a  point  verti¬ 
cally  above  the  target,  some  of  the  computed  values  will  go 
to  infinity,  causing  the  module  to  fail.  The  surface-cutoff 
function  does  not  automatically  turn  itself  off  and  the 


It  has  already  been  mentioned  that  setting  the  IFLU 
switch  equal  to  1  on  the  B-l  record  of  the  STAGS1  run  will 
generate  two  data  files  FOR003.DAT  and  FOR004.DAT  which  are 
specific  to  USA.  These  are  produced  in  addition  to  the 
FOR002.DAT  and  output  files  which  are  normally  generated  by 
a  STAGS1  run . 

Each  of  the  three  data  files  FOR002.DAT,  FOR003.DAT, 
and  FOR004.DAT  must  be  renamed  after  the  STAGS1  run  for  use 
in  USA.  The  first,  FCR^02.DAT,  should  simply  be  renamed  for 
the  file  subdirectory  containing  the  USA  files  (e.g., 

(*. USA) F0R002.DAT).  FOR003.DAT  and  FOR004.DAT  must  be 
renamed  to  correspond  to  the  GRDNAM  and  MASNAM  files  speci¬ 
fied  in  the  FLUMAS  and  AUGMAT  modules,  respectively.  This 
may  be  done  by  entering  a  simple  RENAME  command  on  the  VAX. 
The  user  will  then  be  queried  regarding  the  name  from  which 
and  to  which  they  are  being  changed. 

The  GRDNAM  file  (the  renamed  FOR003.DAT  file)  contains 
the  global  coordinates  of  the  structural  grid  points  gener¬ 
ated  for  the  model  by  STAGS1.  The  FLUMAS  module  then  uses 
this  information  to  construct  a  fluid  mass  matrix  (FLUNAM) 
as  well  as  generate  the  fluid  mesh  data  (GEONAM)  and  trans¬ 
formation  coordinates  that  relate  the  structural  and  fluid 
degrees  of  freedom  on  the  wetted  surface.  The  FLUNAM  and 
GEONAM  data  files  will  then  be  used  in  the  AUGMAT  module  to 
create  a  final  augmented  mass  matrix. 


is  used  to  integrate  across  the  mesh,  STAGS  will  give  erro¬ 
neous  results  (or  none  at  all)  if  that  mesh  is  not  located 
on  the  midplane  of  the  shell  thickness.  ECZ  must  corre¬ 
spondingly  be  set  to  0,  and  IPLAST  set  equal  to  1.  No 
further  inputs  need  be  made  on  this  record  unless  one  of  the 
shell  units  is  located  on  a  plane  of  symmetry.  If  that  is 
the  case,  STAGS  will  not  compute  lateral  deformations 
without  further  information  about  how  the  shell  is  to 
respond.  This  information  is  provided  in  the  form  of  minute 
randomly  generated  displacements  or  "random  imperfections" 
of  the  shell  nodes  from  the  midplane.  To  trigger  input  of 
these  imperfections,  IRAMP  should  be  set  equal  to  1,  and  the 
appropriate  entries  made  on  the  M-6  and  M-7A  records. 

The  user's  manual's  instructions  and  flow  patterns  for 
records  N  through  R  are  clear  and  it  is  believed  that,  with 
possible  references  to  the  input  files  in  Appendix  C,  the  user 
should  have  no  difficulty  completing  the  STAGS  input  deck. 

As  stated  in  the  section  on  organization  of  the  code, 

USA  is  divided  into  four  modules.  Each  of  these  may  be  run 
independently  once  the  needed  data  files  have  been  produced 
by  the  previous  module.  Except  for  the  required  access  to 
those  generated  data  files,  each  module  is  run  entirely  in 
stand-alone  fashion.  Since  modular  input  information  is 
brief  and  error  statements  are  generally  good,  this  permits 
debugging  without  affecting  earlier  inputs  and  generates  a 
high  degree  of  confidence  in  the  modular  output. 


similarly  intended  to  be  a  labor  saving  device.  A  caveat 
should  be  observed,  however,  that  option  1  should  only  be 
used  if  the  side  connecting  two  shells  is  uniquely  defined 
(i.e.,  no  straight  lines).  It  is  strongly  recommended  that 
options  2  or  3,  which  allow  more  user  control  by  specifying 
boundaries  through  specific  global  reference  points, be  used 
in  preference  to  option  1. 

An  option  not  listed  in  the  user's  manual,  but  which 
could  conceivably  be  useful  for  construction  of  cylindrical 
models  is  option  4.  This  allows  the  user  to  specify  the 
location  of  a  unit  origin  in  the  global  coordinate  system 
and  then  define  the  boundary  in  terms  of  one  translation  and 
one  rotation  about  the  reference  point. 

It  should  also  be  noted  that  STAGS  will  permit  introduc¬ 
tion  of  shell  geometries  not  included  in  the  eleven  which 
are  indigenous  to  the  program.  This  is  done  through  a  user 
written  subroutine  which  may  be  specified  by  choosing  ISHELL 
option  1  on  the  M-l  record. 

If  the  ISHELL  =  1  option  is  chosen,  the  shell  geometric 
properties  normally  input  in  record  M-2  will  be  input  in  the 
user-written  subroutine  LAME  as  stated  above.  Otherwise, 
shell  geometries  are  chosen  from  those  provided  by  the  code. 
Property  values  for  these  geometries  should  be  input  as 
specified  in  the  figures  of  Appendix  B.  If  plastic  deforma¬ 
tions  are  expected,  particular  attention  must  be  paid  to 
inputs  to  the  M-5  record.  First,  since  Gaussian  quadrature 


TIMINT  run;  which,  as  has  already  been  mentioned,  is  not 


insignificant  for  the  plastic  case.  The  default  value  for 
plasticity  is  3.  This  may  be  acceptable  for  models  where 
the  membrane  stress  governs,  but  will  not  be  so  if  there  is 
marked  bending  stress. 

Another  change  that  should  be  made  for  runs  including 

plastic  responses  is  to  enter  a  value  of  1  or  2  for  the  LSO 

entry  of  the  K-2  record.  If  this  is  not  done  the  USA-STAGS 

code  will  successfully  complete  an  entire  run,  but  the  post- 

1 8 

processor,  POSTP ,  will  not  print  stress  or  strain  data. 

The  sequence  of  records  M  through  R  defines  the  proper¬ 
ties  and  loads  on  the  shell  units.  As  is  noted  most  explic¬ 
itly  in  the  user's  manual,  these  records  are  read 
sequentially  for  each  unit.  All  data  in  the  sequence  will 
be  read  before  the  M-l  record  for  the  next  unit  is  read. 

The  individual  shell  geometries  and  the  defining  parame¬ 
ters  by  which  they  are  to  be  connected  to  other  shells  (that 
is  the  means  by  which  they  are  to  be  oriented  to  the  global 
geometry)  are  specified  in  record  M-l.  IGLOBE  option  0  is 
useful  if  the  local  coordinates  of  the  unit  already  corre¬ 
spond  precisely  to  the  global  coordinates.  Option  1  is 


18, 


Actually,  this  may  not  be  entirely  true.  In  order  to 
get  any  useful  response  whatsoever  out  of  POSTP,  it  has  been 
necessary  to  "trick"  the  code  into  producing  the  desired 
stresses  and  strains  by  running  the  module  with  an 
F0R002.DAT  file  created  by  a  STAGS1  run  for  the  elastic 
case.  This  has  been  true  whether  or  not  LSO  has  been 
flagged,  with  no  noticable  difference  in  the  results. 
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must  be  consistent  with  the  elastic  modulus  entered  on  the 
1-2  record.  That  is,  the  slope  of  a  line  drawn  from  rhe 
origin  to  that  first  data  entry  should  be  the  elastic 
modulus.  Also,  while  specifying  plastic  behavior 
may  actually  decrease  the  run  time  for  the  TIMINT  module  by 
reducing  the  number  of  iterations  required  for  each  time 
step,  it  can  increase  the  required  computer  space  by  more 
than  100%.^ 

If  plasticity  is  to  be  included  in  the  model,  there  are 
also  some  changes  that  should  be  made  in  the  K  records. 

STAGS  determines  the  plastic  responses  of  the  shell  by 
integrating  through  the  thickness.  The  input  NLIP  in  the 
K-2  record  establishes  the  number  of  integration  points 
through  each  material  layer.  Since  Simpson's  rule  is  used, 
this  number  must  be  odd .  NLIP  may  lie  between  0  and  9 
(inclusive),  but  some  care  must  be  exercised  in  the  choice, 
since  the  larger  numbers  will  add  to  the  length  of  the 


t 


17 

A  normal  run  for  a  model  with  about  300  nodes  seems  to 
be  on  the  order  of  about  two  to  two-and-a-half  hours  of  CPU 
time  for  strictly  elastic  behavior.  If  the  elastic  range  is 
exceeded  but  plastic  behavior  is  not  specified,  these  runs 
can  easily  continue  for  15  to  20  hours  as  the  code  attempts 
repeated  iterations  for  each  time  step  (the  user  can  control 
this  somewhat  by  specifying  the  maximum  allowable  time  for 
any  iteration),  and  the  full  number  of  specified  time  steps 
may  not  be  attainable.  Inclusion  of  plasticity  will  again 
reduce  the  run  time  to  roughly  three  to  four  hours.  Output 
of  stresses  beyond  the  yie.M  point  of  the  material,  however, 
requires  much  larger  solution  data  (FOR022.DAT)  files  for 
both  the  TIMINT  and  POSTP  runs.  Typically,  these  have  been 
on  the  order  of  45,000  and  35,000  blocks  respectively  for 
the  pressure  levels  used  here. 
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originally  done.  In  an  attempt  to  make  the  USA-STAGS 

output  correspond  directly  to  the  strain  gage  outputs  from 

underwater  testing,  the  plates  were  modelled  with  the  meshes 

placed  on  the  surfaces  of  the  plates  which  corresponded  to 

the  surface  of  strain  gage  mounting.  The  half-thickness  of 

the  plate  was  then  input  as  the  ECZ  value  on  the  M-5  record, 

thereby  defining  the  true  mid-plane  for  the  code.  This 

procedure  worked  extremely  well  for  elastic  runs.  As  stated 

in  the  previous  section,  however,  the  Gaussian  quadrature 

that  STAGS  uses  to  calculate  the  plastic  responses  across 

the  mesh  surface  is  incapable  of  handling  the  integrations 

25 

at  any  but  the  mid-plane. 

For  each  of  the  test  plates,  maximum  advantage  has  been 
taken  of  geometric  symmetries  in  order  to  reduce  the  size  of 
the  finite  element  model.  For  the  tee-stiffened  plate  this 


25 

This  important  point  is  reiterated  here  because  it  may 
not  be  immediately  apparent  from  the  STAGS  run.  The  exper- 
ence  gained  here  seems  to  indicate  that  when  a  negative 
eccentricity  is  input  to  the  code,  perfectly  reasonable  (but 
wrong!)  strain  approximations  are  obtained.  The  first  indi¬ 
cation  of  error  in  this  method  for  correlating  STAGS 
and  strain  gage  data  arose  when  a  run  with  a  positive 
eccentricity  was  attempted.  Under  those  conditions,  STAGS 
produced  no  strain  data  at  all.  This  inability  of  STAGS  is 
not  incapacitating,  but  it  is  inconvenient  for  it  means  that 
the  data  generated  in  the  computer  run  will  have  to  be 
converted  to  correspond  to  the  strain  gage  information.  In 
theory  this  is  simply  a  matter  of  multiplying  the  plate 
half-thickness  by  the  curvature  calculated  in  the  run  and 
then  adding  or  subtracting  that  from  the  mid-plane  strain. 
Given  the  massive  size  of  the  STAGS  output,  however,  this 
means  that  either  special  processors  will  have  to  be 
employed  or  a  conversion  program  will  have  to  be  written. 
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presents  no  further  complications  since  the  resulting 
quarter  model  is  non-symmetr ic .  For  the  plate  with  the 
longitudinally  oriented  rectangular  stiffener,  however,  the 
result  of  half  plate  modelling  is  that  the  stiffener  itself 
lies  on  a  plane  of  symmetry.  This  presents  a  problem  for 
STAGS.  Although  deformations  along  the  symmetric  plane  will 
be  calculated,  no  lateral  displacements  will  be  recognized 
since,  as  far  as  STAGS  is  concerned,  forces  on  either  side 
are  equal.  This  both  affects  the  validity  of  the  displace¬ 
ments  that  are  calculated  and  causes  a  failure  to  indicate 
stiffener  tripping.  STAGS  solves  this  problem  by  allowing 
the  user  to  either  input  initial  geometric  imperfections  in 
the  shell  reference  surface  or  specify  that  the  code 
generate  random  imperfections.  As  stated  in  the  previous 
section,  this  is  done  by  flagging  the  IWIMP  and  I RAMP  vari¬ 
ables  on  the  M-5  record  and  then  following  the  flow  direc¬ 
tions  listed  in  the  manual  through  the  M-7A  record.  For  the 

rectangularly  stiffened  plate  used  here,  imperfections  in 

2  0 

the  first  three  harmonics  of  0.001  inch  were  specified. 

Finally,  although  not  strictly  related  to  plate  discre¬ 
tization,  there  is  one  other  point  not  mentioned  in  the 
STAGS  user's  manual,  but  which  is  essential  if  the  computer 
code  is  to  successfully  compute  plastic  def ormat ions .  If 

26 

See  the  STAGS1  inputs  for  the  rectangular  stiffened 
plate  in  Appendix  C. 


plasticity  is  to  be  included,  the  variable  IPLST  on  the  1-1 

record  must  be  set  to  some  value  between  0  and  4  to  indicate 

01 

which  of  four  plastic  strain  theories  should  be  used.  A 
non-zero  positive  value  must  also  be  specified  for  the  vari¬ 
able  NESP .  This  NESP  value  specifies  the  number  of  data 
points  on  a  stress-strain  curve  for  the  plate  material  which 
will  be  input  to  the  1-3  record.  There  are  two  points  that 
the  user's  manual  fails  to  make  clear  about  these  stress- 
strain  points.  First,  this  curve  requires  stress  and  engi¬ 
neering  strains  beyond  the  yield  point.  Second,  the  first 

point  of  the  curve  must  correspond  closely  to  the  modulus  of 

2  8 

elasticity  of  the  material  specified  in  the  1-2  record. 

All  the  material  and  stress-strain  characteristics  for  the 
6061-T6  aluminum  used  to  manufacture  the  test  plates  here 
were  taken  from  the  Military  Standardization  Handbook, 
Metallic  Materials  and  Elements  for  Aerospace  Vehicle 
Structures  [Ref.  14:  p.  3-184]. 

27 

The  default  value,  0,  automatically  triggers  use  of 
the  White-Bessel ing  theory,  as  does  the  set  value  1. 

28 

That  is,  the  slope  from  the  origin  to  the  first  point 
specified  must  be  the  modulus  of  elasticity.  The  White- 
Besseling  theory  will  apply  approximations  to  attempt  a 
patch-up.  If  the  discrepancy  is  large,  however,  continuity 
between  the  elastic  and  plastic  responses  will  be  lost  and 
the  resulting  output  will  be  in  error. 


V.  TESTING  THE  CODE 


A.  MODAL  TESTING 

Determination  of  a  structure's  mode  shapes  and  natural 
frequencies  can  provide  a  first  test  of  a  finite  element 
code's  ability  to  predict  the  dynamic  responses  of  an 
excited  system.  Perhaps  of  even  more  importance  to  the 
purposes  of  this  study,  it  can  be  an  excellent  check  of  the 
structural  data  input  to  the  code.  Failure  to  provide 
proper  connectivity  between  shell  units  or  inaccuracies  in 
the  input  dimensions  or  material  properties  will  show  as 
gross  discrepancies  in  the  modal  analysis. 

The  Naval  Postgraduate  School  is  fortunate  to  be  one  of 
the  few  institutions  to  have  a  Hewlett-Packard  System  5451C 
Fourier  Analyzer  with  modal  testing  capability. 

The  technique  used  by  the  H-P  5451C  is  in  many  ways 
similar  to  the  solution  processes  used  by  the  finite  element 
codes  themselves  for  it  involves  solving  basically  the 
same  governing  equation  of  motion  [Ref.  15]. 

Mx  +  Cy  +  Ky  +  f  =  0  (eqn  5.1) 

Since  both  x  and  f  are  functions  of  time,  the  transfer 
function,  h(x)»  may  be  found  by  taking  the  Laplace  trans¬ 
form  of  Equation  5.1 


B(s )X(s )  =  F(s)  ( eqn  5.2) 

where , 

B(s)  =  Ms^  +  Cs  +  K  (eqn  5.3) 

and , 

H(s)  =  B(s)-1  (eqn  5.4) 

If  the  exciting  forces  and  displacements  are  known,  H(s) 
may  be  readily  found  by  solving  the  equation  which  results 
from  substituting  Equation  5.4  into  Equation  5.2,  namely 

H(s)F(s)  =  X(s)  (eqn  5.5) 

For  an  oscillating  system  of  order  n,  the  transfer  func¬ 
tion  will  have  2n  poles  occurring  in  complex  conjugate  pairs 
Each  pair  of  poles  will  cause  a  mode  of  vibration  in  the 
structure.  Each  of  these  poles  is  a  complex  number  which 
may  be  expressed  as 

=  -a^  +(^i  (eqn  5.6) 

where  is  the  damping  coefficient  and  the  natural 
frequency  for  that  mode. 

Once  and  have  been  determined,  the  resonant 
frequency  (in  radians  per  second)  and  the  damping  ratio 


C  follow  directly. 


ak  + 


(eqn  o .  7 ) 


(eqn  5.8) 


The  exciting  forces  and  displacements  that  must  be 

provided  to  the  system  may  be  found  in  a  number  of  ways. 

The  method  chosen  here  was  to  excite  the  plate  using  a 
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modally  tuned  impulse  hammer.  An  accelerometer  was 
attached  to  the  plate  to  record  the  resulting  accelerations, 
from  which  displacements  were  found  by  numerical 
integration . 

The  two  fundamental  assumptions  of  this  test  procedure 
[Ref .  15 ]  are  that : 

1.  Modal  frequency  and  damping  are  constants  for  any 
transfer  function  taken  from  the  structure. 

2.  Modes  of  vibration  can  be  excited  from  anywhere  on 

an  elastic  structure  except  at  their  node  lines  (where 
no  excitation  is  possible). 

It  should  therefore  be  possible  to  describe  all  modes  of 
interest  by  simply  choosing  a  sufficient  number  of  excita¬ 
tion  points  on  the  structure  and  by  placing  these  points 
close  enough  together  to  minimize  the  possibility  of 
replicat ion . 

For  the  two  plates  used  in  these  tests  (the  rcctangular- 
and  tee-stiffened  test  models),  nine  points  of  excitation 


PCB  model  086B03. 


PCB  model  302A07. 


were  chosen  on  the  back  of  each  plate.  Care  was  taken  to 

32 

avoid  symmetric  placements  of  the  points.  The  accelerometer 
was  fixed  at  a  single  location  arbitrarily  chosen  but  spaced 
away  from  the  clamped  boundaries. 

When  any  location  was  tapped  with  the  modal  hammer,  the 
resulting  forces  were  fed  back  through  the  hammer  to  the 
Fourier  Analyzer;  as  were  the  accelerations  from  the  accel¬ 
erometer.  The  Fourier  Analyzer  then  computed  the  displace¬ 
ments  and  solved  the  Laplace  transform  to  create  the  plot 
shown  in  Figure  5.1.  Each  of  the  peaks  in  this  representa¬ 
tion  corresponds  to  the  poles  of  the  transfer  function  for  a 
single  point  of  excitation.  By  then  combining  the  plots  for 
all  points,  an  accurate  representation  of  the  system  modes 
was  obtained.  It  was  then  simply  a  matter  of  analyzing  each 
of  the  poles  in  descending  order  of  magnitude  to  whatever 
total  number  of  modes  is  of  interest  to  the  user  to  deter¬ 
mine  the  dominant  modes  of  the  system. 

B.  UNDERWATER  SHOCK  TESTING 

Underwater  explosion  testing  for  all  studies  conducted 
in  this  series  to  date  has  been  done  at  the  West  Coast  Shock 
Test  Facility  (WCSF),  located  at  the  old  Hunter's  Point 
Naval  Shipyard  near  San  Francisco,  California.  This  is  a 

31 

The  plates  were  placed  face  downward  upon  a  solid 
marble  slab. 
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The  actual  number  and  spacing  of  the  excitation  points 
seems  to  be  primarily  a  matter  of  experience. 

88 
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Navy  activity  reporting  to  the  Supervisor  for  Shipbuilding, 
Conversion,  and  Repair,  San  Francisco.  As  the  only  licensed 
facility  of  its  type  on  the  West  Coast,  they  have  gained  a 
considerable  amount  of  experience  in  conducting  underwater 
explosion  testing;  although  most  of  it  has  been  large  scale 
work  in  support  of  the  Navy's  shock  qualification  program 
rather  than  experimentation. 

1 .  The  Test  Platform 

In  order  to  provide  the  air  backed  condition  that 
was  required  to  simulate  a  ship's  hull  and  to  ensure  rigidly 
clamped  boundary  conditions,  the  test  plates  used  in  this 
series  of  studies  we  e  securely  bolted  to  a  heavy  steel 
backing  structure  (Fig.  5.2).  This  structure  consisted  of  a 
box  constructed  of  li  inch  thick  A6  structural  steel.  A 
wide  flange  was  welded  around  the  open  side  of  the  box  to 
which  the  edge  of  the  plate  was  bolted.  An  O-ring  gasket 
was  fit  into  a  channel  machined  into  the  surface  of  the 
flange  to  provide  watertight  integrity.  Access  to  instru¬ 
mentation  on  the  inside  of  the  box  was  allowed  by  a  subma- 
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rine  style  penetrator  affixed  through  the  bottom. 

To  perform  the  actual  testing,  the  plate  and  backing 
structure  were  suspended  by  heavy  steel  cables  below  two 
large  pneumatic  floats  in  the  configuration  developed  in 


For  a  more  complete  description  of  this  backing  struc¬ 
ture  and  the  rationale  used  in  designing  it,  see  Reference 

3,  pp.  88-94. 


Chapter  2  (Fig.  5.3).  The  cables  were  cut  to  length  to 
assure  that  the  nearest  edges  of  the  test  panel  stiffeners 
were  maintained  at  depths  as  close  to  thirteen  feet  below 
the  surface  as  possible.  Once  the  entire  apparatus  had  been 
towed  into  the  bay,  the  explosive  charge  was  lowered  into 
position  nine  feet  above  the  stiffeners  in  line  vertically 
with  the  center  of  the  plate. 

2 .  Instrumentat ion 

In  his  work,  which  has  already  been  cited  several 
times,  LT  Rentz  outlined  the  basic  theory  governing  the  test 
of  the  flat  plate  and  began  developing  the  procedures  which 
will  be  enlarged  upon  here. 

Essentially,  the  information  sought  for  comparison 
with  the  computer  simulation  codes  being  tested  has  been: 

1.  Fluid  pressures  at  the  plate  and  in  the  freefield. 

2.  Strain  in  the  plate. 

3.  Final  deflections  in  the  plate. 

Instrumentation  used  to  gather  this  data  was  as  follows: 

a.  Pressure  Measurements 

Three  pressure  gages  were  used  in  the  tests 
conducted  here.  Two  of  these  were  located  just  as  they  had 
been  in  the  previous  study  [Ref.  3:  p.  46].  One  was  located 
one  foot  above  the  center  of  the  plate  and  the  other  was 
clamped  to  the  aluminum  block  bolted  to  the  side  of  the 
plate.  The  forward  edge  of  this  second  pressure  gage  was 


positioned  so  as  to  be  on  an  arc  of  the  same  radius  as  the 

standoff  distance'. 


The  pressure  recordings  that  had  been  obtained 
from  these  two  gages  had  been  so  subject  to  interference 
from  the  plate  structure  and  the  rigging  of  the  test  plat¬ 
form  that  the  futility 'of  conducting  further  tests  without  a 
true  free  field  pressure  gage  was  readily  obvious.  Accor¬ 
dingly,  although  the  first  two  gages  were  retained  as 
backups,  a  third  gage  was  positioned  on  an  arc  equal  to  the 
standoff  distance  but  removed  five  feet  from  the  centerline 
of  the  plate  (Fig.  5.4).  As  will  be  discussed  in  the 
following  chapter,  the  results  obtained  from  this  gage 
during  the  test  shot  on  the  tee-stiffened  plate  were 
excellent . 

All  tests  in  this  series  utilized  one-quarter 

inch  tourmaline  crystal  pressure  gages  rated  to  10,000  psi. 

34 

The  response  ratio,  R  ,  is  a  function  of  the  decay 

P 

constant  for  the  charge,  9,  divided  by  the  transit  time  of 
the  shock  wave  across  the  gage,  t^  [Ref.  16].  For  the 
eight  pound  TNT  charge  chosen  as  the  standard  for  these 
studies , 

9  =  0.15  msec 

34 

Defined  as  the  ratio  of  the  apparent  peak  pressure  to 
the  actual  peak  pressure. 
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For  a  quarter  inch  gage, 


= 


0.25  in 


'D 


60,000 


m 


sec 


-6 


=  4.167  x  10  v  sec 


Applying  the  correction  factor  for  the  oil-filled  protective 
boot  which  surrounds  the  crystal, 


,-6 


1.5  tD  =  6.25  x  10  sec 


Thus  , 


1.5  t. 


=  24 


From  the  chart  on  page  17  of  Reference  16. 


R  =  0.9733 
P 

This  indicates  that  a  very  good  correlation  can  be  expected 
between  the  actual  and  measured  pressures,  and  is  also  the 
best  indication  that  the  correct  gage  size  is  being  used, 
b.  Strain  Measurements 

Perhaps  the  single  most  critical  aspect  of 
obtaining  successful  data  from  underwater  explosion  testing 
of  the  type  conducted  here  is  the  choice  and  application  of 
the  strain  gages. 


For  the  tee-stiffened  plate,  twenty  SR-4  type 
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FAE-25-12-S13ET  strain  gages  were  attached  as  shown  in 
Figure  5.5.  These  were  the  same  type  used  in  Reference  3; 
and,  although  rated  to  only  2000  microstrain,  it  had 
been  demonstrated  that  their  range  could  be  extended  suffi¬ 
ciently  to  provide  useable  data  over  the  very  earliest 
increments  of  the  shock  response.  As  in  the  earlier  test, 
BLH's  SR-4  Permabond  910  was  used  as  the  fixative,  and  the 
attached  strain  gages  and  their  wiring  was  covered  with 

PR-1422  A-l  Base  Compound,  manufactured  by  Products  Research 
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and  Chemical  Corporation  of  Glouster  City,  NJ .  This  base 

compound  was  also  used  on  the  plate  with  the  rectangular 

stiffener  where  19  MX  type  PAHE-250BA-350EN  strain  gages37 

(Fig.  5.6)  rated  to  approximately  200,000  microstrain  were 

used  in  an  attempt  to  increase  the  time  span  over  which 

information  could  be  gathered.  Also,  at  the  recommendation 

of  the  distributor,  BLH  type  EPY-150  epoxy  adhesive  was 

38 

used  instead  of  the  Bondo  910. 


Manufactured  by  BLH  Electronics  of  Waltham,  MA . 

30 

This  was  done  to  provide  and  smooth  compliant 
covering  that  would  keep  the  gages  and  wiring  from  being 
torn  off  with  the  passage  of  the  shock  wave,  and  to  seal  the 
gages  from  saltwater  once  the  rig  was  submerged. 

37 

Manufactured  by  Experimental  Stress  Analysis  Products. 
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In  theory,  the  epoxy  provides  a  stronger  bond  than 
does  the  Bondo,  which  is  a  contact-type  adhesive. 
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In  practice,  although  some  useful  data  was 
obtained  from  each  test,  neither  was  a  particularly 
successful  source  of  strain  information.  The  number  of 
gages  lost  far  exceeded  the  number  that  performed. 

In  recapitulation,  a  number  of  faults  were 
found  for  which  the  following  corrective  measures  have  been 
proposed. 

1.  Although  the  MX  type  strain  gages  greatly  exceed 
the  strain  requirements  expected  in  these  tests, 
their  use  is  strongly  recommended.  Little  or  no 
loss  of  accuracy  is  experienced  in  the  lower  strain 
regions  through  use  of  the  high-elongation  gages, 
whereas  mid-range  gages  tend  to  peak  right  at  the 
strain  ranges  expected. 

2.  On  any  strain  gage,  the  weakest  points  are  the  sold¬ 
ering  tabs  to  which  the  gage  wires  are  attached. 

These  are  brittle  and  will  break  if  any  bending 
stress  is  applied.  All  motion  should  be  restricted 
to  the  plane  of  the  gage.  This  means  that  the  strain 
gage  wires  must  be  extremely  flexible  so  that  any 
bending  motion,  either  during  installation  or  as  a 
result  of  plate  deflections,  is  not  transmitted  to 

the  gage  tabs.  Rated  strain  gage  wire’  is 
39 

recommended . 
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On  the  tee-stiffened  plate,  28  gage  wire  made  specifi¬ 
cally  for  strain  gages  was  used.  This  proved  to  be  quite 
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3.  Even  with  flexible  strain  gage  wire,  a  determined 


force  along  the  length  of  the  wire,  as  would  be  experi 

enced  during  plate  deformation  will  break  the  tab, 

the  wire,  or  the  solder  connection  between  the  two. 

Although  a  plastic  medium,  the  PRC  Base  Compound 

used  to  protect  the  gages  and  wires  was  not  flexible 

40 

enough  to  allow  the  wires  to  expand  with  the  plate. 

A  suggested  solution  to  this  problem  is  to  lay  the 
gage  wires  with  some  s-curves  to  allow  elongation  as 
the  plate  deforms.  The  wires  should  not,  of  course, 
be  otherwise  affixed.  A  piping  of  Barrier  WD ,  a 
tacky,  semi-fluid  substance  manufactured  by  BLH  which 


satisfactory  as  far  as  the  gages  themselves  were  concerned, 
but  presented  definite  problems  on  the  human  side  of  the 
equation.  The  wire  was  so  fine  that  it  was  impossible  to 
work  with  it.  It  also  proved  to  be  highly  vulnerable  to 
damage  while  the  test  structure  was  being  transported  or 
handled.  For  the  following  shot,  the  attempt  was  made  to 
use  20  gage  bell  wire.  This  withstood  handling  magnifi¬ 
cently,  but  was  too  stiff  for  the  soldering  tabs.  Annealed 
22  or  24  gage  wire  designed  specifically  for  use  with  strain 
gages  appears  to  be  the  best  possible  choice. 
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This  base  compound  seems  to  be  an  extremely  protective 
medium.  A  couple  of  cautions  should  be  issued  governing  its 
use,  however.  First,  while  the  compound  must  be  used  within 
one  hour  of  mixing  base  and  accelerator,  enough  of  that  time 
should  be  dedicated  to  stirring  to  ensure  that  the  compound 
is  thoroughly  mixed.  Although  this  observation  is  not 
completely  confirmed,  it  appears  that  insufficient  mixing 
leaves  a  residual  which  can  permeate  wiring  insulation  and 
actually  destroy  very  fine  copper  wires.  Second,  the 
compound  is  not  compatible  with  the  Gage  Cote  lacquer 
normally  painted  over  strain  gages  as  a  sealant.  If  the 
base  compound  is  to  be  applied  directly  to  the  strain  gages, 
Gage  Cote  should  not  be  used. 


looks  somewhat  like  a  combination  of  beeswax  and 

petroleum  jells-,  could  then  be  applied  over  both  gages 
41 

and  wires.  This  piping  of  Barrier  WD  should  then  be 

covered  by  a  layer  of  aluminum  foil  and  then  the  whole 

overlayed  by  the  usual  PBC  Base  Compound.  The  net 

effect  is  an  armored  "tube"  of  Barrier  WD  through  which 

the  wires  can  run  when  the  plate  deforms. 

4.  The  fixative  used  to  attach  the  strain  gages  should, 

itself,  be  a  subject  for  careful  consideration.  Even 

though  the  epoxy  is  supposed  to  be  the  stronger  medium, 

experience  gained  here  indicates  that  it  is  not  as 

satisfactory  as  the  Bondo  for  explosion  testing.  The 

postulated  reason  for  this  seeming  discrepancy  is  that 

the  epoxy  places  a  thicker  layer  of  adhesive  between 

gage  and  plate  material.  Given  the  extremely  rapid 

loadings  experienced,  intense  shears  are  created  in  the 

thicker  substance  and  the  epoxy  literally  cracks  and 

spalls  off.  The  Bondo  910  therefore  seems  to  be  the 

42 

preferable  adhesive. 
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A  bakers  pastry  tube  works  very  nicely  to  do  this. 
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No  failures  due  to  loading  were  noted  in  the  strain 
gages  attached  with  Bondo.  Therefore,  Even  if  the  cracking 

problem  had  not  occurred  in  the  epoxy,  the  greater  separation 
between  gage  and  plate  material  that  the  epoxy  creates  would 
argue  in  favor  of  the  Bondo. 


Plate  Deflections 


c . 

Plate  deflections  were  by  fa,’  the  most  reliable 

source  of  data  that  could  be  obtained  as  a  result  of  testing 

since  no  electronics  were  involved.  To  get  these,  precise 

measurements  were  made  at  specified  points  on  the  plate  both 

before  and  after  the  test  shot.  A  jig  was  devised  which 

took  advantage  of  the  rigid  plate  boundaries  to  provide  a 

reference  table  for  a  dial  indicator  depth  gage  (Fig.  5.7). 

All  measurements  were  made  while  the  test  plate  was  bolted 

onto  the  backing  structure  to  ensure  that  plate  warpage  did 

43 

not  affect  the  readings. 

C.  DATE  RETRIEVAL 

Signals  from  the  gages  mounted  on  and  about  the  test 
plate  were  fed  back  through  shielded  cable  to  a  common 
floating  terminal  box.  From  there,  they  were  passed  to  an 
instrumentation  shack  located  on  one  of  the  jetties  defining 
the  small  bay  in  which  the  testing  was  conducted.  A  sche¬ 
matic  of  the  setup  is  shown  in  Figure  5.8.  Calibration 
levels  were  applied  to  all  pressure  signals  by  Endevco 
Model  4470  Signal  Conditioners.  These  same  signal  condi¬ 
tioners  were  used  to  provide  calibration  for  all  but  twelve 
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Deflections  have  already  been  given  for  the  rectan¬ 
gular  stiffened  plate  in  Figures  3.7  and  3.8.  No  post-shot 
measurements  were  made  of  the  tee-stiffened  plate  because  of 
the  heavy  damage  which  was  inflicted  when  the  plate  was  hit 
by  one  of  the  supporting  cables. 
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of  the  strain  gage  channels.  Those  remaining  twelve 
channels  were  calibrated  by  two  six-channel  Honeywell  Model 
B2-6  Universal  Bridge  Balance  units.  The  strain  gage  signal 
were  then  amplified  by  Bell  and  Howell  CEC  1-168  DC  ampli¬ 
fiers  before  final  recording.  All  signals  were  recorded 
on  two  twelve  channel  Ampex  FR  1300  tape  recorders. 

Calibration  level's  for  the  tourmaline  crystal  pressure 
gages  were  all  set  at  the  10,000  psi  theoretical  maximum. 

For  the  strain  gages,  however,  attempts  were  made  to  use  the 
computer  codes  to  provide  estimates  of  expected  strains 
using  the  pressure  and  decay  constant  calculated  from 
Equations  2.2  and  2.3.  This  was  especially  critical  for 
the  strain  levels  for  the  tee-stiffened  plate  where  some  of 
the  strain  gages  were  expected  to  be  operating  at  or  beyond 
their  rated  limits. 

Unfortunately,  difficulties  encountered  in  learning 

to  use  the  USA-STAGS  code  prevented  runs  from  being  made 
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before  the  shot  dates. 

Advantage  could  have  been  taken  of  earlier  estimates 
used  for  the  plate  tested  in  Reference  3  if  the  charge 
weights  had  remained  the  same;  however,  such  was  not  the 
case . 

Since  8  pound  TNT  charges  are  listed  as  Navy  standard 
stock  items  and  order  had  been  placed  well  in  advance,  no 
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In  fact  no  successful  runs  were  made  until  long  after 
all  experimentation  had  been  completed. 


problems  were  anticipated  in  obtaining  charges  of  this  size. 
As  the  date  of  the  first  test  neared,  however,  it  became 
obvious  that  no  8  pound  charges  would  be  available.  The 
nearest  alternate  that  could  be  provided  in  time  for  the 
scheduled  test  was  a  9i  pound  TNT  charge. 

It  can  perhaps  be  appreciated  that  dimensional  changes 
to  the  test  platform  are  not  simply  a  matter  of  moving  the 
charge  location.  As  charge  weights  are  changed,  the  stand¬ 
off  distance  and  the  charge  depth  both  must  be  changed; 

which  involves  cutting  new  supporting  cables  and  recali-, 
brating  flotation  depths.  All  of  which  is  extremely  time 
consuming.  When  the  personnel  at  WCSF  were  contacted  about 
the  possibility  of  changing  the  configuration  for  a  9A 
pound  charge,  they  advised  that  it  co-uld  not  be  done  in  time 
for  the  scheduled  tests. 

The  choice  for  the  first  test  was  clear,  either  run  the 
test  w.  th  the  existing  platform  and  a  9i  pound  charge,  or 
abort  the  experiment.  It  was  determined  that  more  benefit 
could  be  gained  from  conducting  the  test. 

The  values  for  Pq  and  9  found  from  Equations  2.2  and 
2.3  for  a  9£  pound  charge  are: 

Pq  =  4081.7  psi,  9  =  0.161  msec 

This  initial  pressure  was  markedly  higher  than  that  for 
an  8  pound  charge.  The  maximum  bubble  radius,  too,  was 


Before  the  test  platform  was  placed  in  the  water,  all  o 
the  strain  gage  signals  were  balanced  with  their  calibratio 
resistors,  setting  their  signals  to  zero.  As  soon  as  the 
platform  was  placed  in  the  water,  the  change  in  temperature 
effected -a  corresponding  change  in  the  resistance  of  the 
strain  gages.  This  in  turn  changed  the  balance  of  the 
bridge.  This  was  an  ongoing  problem  throughout  the  period 
of  the  tests;  many  of  the  strain  gages  never  were  brought 
into  proper  balance.  The  problem  was  exacerbated  by  the 
fact  that  the  strain  gage  and  its  associated  wiring  formed 
an  extremely  long  leg  of  the  balance  bridge  which  was 
subject  to  the  environment  of  the  bay.  The  balancing  leg 
consisted  of  a  calibration  resistor  located  in  the  instru¬ 
ment  shack;  a  very  different  environment  from  that  outdoors 
The  second  of  these  three  problems  is  one  which  is 
readily  correctable  by  the  experimenter  if  he  is  aware  of 
it.  When  the  test  plate  and  backing  structure  were  sent  to 
the  West  Coast  Shock  Test  Facility,  all  strain  gage  leads 
were  connected  to  junctions  in  the  terminal  box  located  on 
the  side  of  the  backing  structure.  Connecting  these  leads 
through  long  cable  runs  was  one  of  the  many  tasks  that  was 
performed  by  WCSF  personnel  between  the  time  that  the  plate 
was  delivered  and  the  date  of  the  test  shot.  Their  experi¬ 
ence  was  primarily  with  the  Navy's  shock  test  program, 
wherein  the  magn  1 1 udos  of  the  strain  outputs  were  of  more 
value  than  the  sign.  Therefore,  little  attention  was  paid 


The  voltage  used  to  detonate  the  explosive  charge  is 
quite  large  (800  volts  DC)  and  sets  up  a  noticeable  induc¬ 
tance  on  all  electrical  cabling  in  the  vicinity.  In  fact, 
on  the  strain  gages,  which  are  very  low  impedance  devices 
and  therefore  relatively  insusceptible  to  ind'Ctive  noise 

this  was  used, as  a  convenient,  if  not  very  accurate,  zero 

*  .  .48 

time  mark. 

The  pressure  gages,  on  the  other  hand,  are  extremely 
high  capacitance,  high  impedance  gages  and  correspondingly 
are  very  much  affected  by  induced  currents.  It  is  therefore 
believed  that  the  close  proximity  of  the  detonation  cable 
was  the  sole  cause  of  the  saturation  observed  on  the  pres¬ 
sure  gages  during  the  second  test.  Be  that  as  it  may,  no 
pressure  history  was  obtained  for  the  rect angular-st i f f ened 
plate,  making  it  necessary  to  use  the  peak  pressure  and 
decay  constant  calculated  from  Equations  2.2  and  2.3  in 
the  subsequent  USA-STAGS  run. 

The  lack  of  a  zero  time  signal  was  not  the  only  defi¬ 
ciency  noted  on  the  strain  gage  histories.  Three  others 
should  also  be  mentioned  by  way  of  warning  for  those 
conducting  similar  studies. 
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This  was  a  major  deficiency  in  the  Way  testing  was 
conducted  at  WCSF.  No  accurate  zero  time  was  provided  on 
any  channel.  The  detonation  signal  was  used  as  a  best 
alternative,  but  since  that  typically  covered  a  span  of 
about  one  millisecond  it  was  almost  useless  when  only  the 
first  few  milliseconds  were  being  studied. 
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the  peak  pressure  for  a  £  pound  TNT  charge  is  far  less 
than  that . 


The  methods  by  which  the  two  tests  were  conducted 
suggested  a  probable  cause  for  the  pressure  gage  saturation. 

Once  all  wiring  was  connected  and  the  gages  tested,  the 
test  platform  was  lowered  into  the  waters  of  the  bay.  Heavy 
halyards  were  attached  to  either  end  of  the  floats,  one  to 
be  used  as  an  outhaul  and  the  other  as  an  inhaul .  Using 
these,  the  platform  was  positioned  in  the  bay  a  safe 
distance  from  the  jetties  before  the  explosive  charge  was 
installed  and  the  test  conducted. 

To  conduct  the  test  with  the  pound  charge,  the 
platform  was  centered  in  the  bay  to  keep  it  as  far  away  from 
the  jetties  as  possible.  There  was  a  heavy  wind  blowing  and 
when  the  cable  used  to  carry  the  firing  current  to  the 
explosive  charge  was  taken  out  to  the  platform  it  was  blown 
well  away  from  the  platform  and  the  instrument  cables. 

Since  it  was  separated  from  the  other  instruments,  the  free- 
field  pressure  gage  cable  was  farthest  of  all  from  the 
detonating  cable. 

When  the  test  shot  was  conducted  on  the  rectangular- 
stiffened  plate,  the  charge  was  sufficiently  small  that  the 
platform  was  kept  much  nearer  pierside  to  make  it  easier  to 


handle.  All  of  the  cabling  was  run  out  to  the  platform 
along  the  inhaul . 
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(3125  psi)(.  36788)  =  1149.6  psi 

which  occurs  at  about  1.56  milliseconds. 

These  values  were  input  to  the  USA-STAGS  code  to  conduct 
a  post-shot  analysis  for  comparison  to  test  data. 

The  peak  pressure  and  decay  constant  obtained  from  the 
test  results  are  much  smaller  than  those  predicted  by 
Equations  2.2  and  2.3.  In  retrospect,  this  is  not  parti¬ 
cularly  surprising.  First,  it  will  be  recalled  that  Equations 
2.2  and  2.3  are  derived  from  empirical  results  and  are  there¬ 
fore  subject  to  wide  variations.  Even  more  pertinent,  it 
was  found  in  Chapter  V  that  the  ratio  of  the  charge  depth  to 
the  maximum  bubble  radius  was  much  less  than  the  .50 
threshold  for  pressure  loss.  The  decreased  pressure  and 
more  gradual  decay  are  exactly  what  would  be  expected  from  a 
charge  that  had  vented  a  significant  portion  of  its  energy 
in  the  first  moments  after  detonation. 

Given  the  good  results  of  the  first  shot,  a  high  degree 
of  optimism  was  established  for  the  results  of  the  pressure 
data  from  the  second,  the  £  pound  charge  on  the  plate  with 
the  longitudinal  rectangular  stiffener.  That  optimism  was 
not  justified. 

Upon  playback  of  the  data  for  the  second  shot  ,  it  became 
clear  that  the  amplifiers  from  all  three  of  the  pressure 
gages  had  become  totally  saturated.  This  was  surprising 
since  the  pressure  gages  were  calibrated  for  10,000  psi  and 
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pressure  and  strain  gages  used  for  the  tests  along  with 
their  associated  calibration  and  time  signals  were  recorded 
by  two  Ampex  FR  1300  tape  recorders.  After  each  of  the  test 
shots,  the  data  tapes  and  one  of  the  tape  recorders  were 
transported  to  the  Naval  Postgraduate  School  where  the  tapes 
were  played-back,  filtered,  and  recorded  on  disk  on  the  H-P 
5451C  Fourier  Analyzer  for  later  analysis  and  plotting. 

As  thorough  a  pretest  analysis  as  was  possible  without 
actual  disassembly  was  conducted  upon  the  equipment  at  the 
West  Coast  Shock  Test  Facility  before  each  shot.  Calibra¬ 
tion  levels  for  all  sensing  devices  were  also  checked 
using  technicians  and  equipment  provided  by  the  Naval 
Postgraduate  School.  Despite  these  precautions,  the  quality 
of  data  obtained  from  the  tests  was  highly  problematical. 

An  excellent  free-field  pressure  history  was  obtained 
from  the  test  shot  conducted  on  the  tee-stiffened  plate 
(Fig.  6.1);  fully  compensating  the  difficulties  encountered 
in  mounting  the  third  pressure  gage  in  a  position  where  it 
would  not  be  affected  by  interference.  As  can  readily  be 
observed,  this  curve  displays  all  of  the  classical  exponen¬ 
tial  decay  characteristics  predicted  in  Chapter  II.  Since 
this  is  so,  the  decay  constant  will  correspond  approximately 
to  the  time  it  takes  for  the  pressure  to  drop  to  ^  (about 
.36788)  of  its  peak  value.  From  Figure  6.1,  the  peak 
incident  pressure  was  3125  psi.  Therefore, 


milled,  nor  was  it  possible  to  ensure  an  absolutely  flush 
mount  between  plate  and  slab. 

When  the  tee-stiffened  plate  was  tested,  the  bolt  holes 

had  already  been  drilled  in  the  frame.  These  bolt  holes 

seemed  to  cause  reflections  which  were  picked-up  by  the 

Fourier  Analyzer  as  additional  modes.  These  may  be  observed 

to  some  degree  in  Table  IV  in  the  repeated  tendency  for  a 

group  of  frequencies  to  cluster  about  a  single  value.  The 

effect  was  far  more  noticeable  on  the  oscilloscope  traces 

obtained  as  each  of  the  excitation  points  on  the  plate  was 
47 

loaded.  When  all  of  the  variables  are  considered,  it  is 
believed  that  the  actual  correspondence  between  the  modal 
test  results  and  the  STAGS2  prediction  is  quite  good  and 
provides  adequate  confirmation  of  both  the  code’s  ability 
to  predict  the  natural  frequencies  of  a  simple  structure 
and  of  the  finite  element  models  input  to  the  code  for  this 
study . 

B.  COMPARISONS  WITH  UNDERWATER  TEST  RESULTS 

As  has  already  been  stated  in  Chapter  V,  all  underwater 
shock  testing  for  this  series  of  studies  was  conducted  at 
the  West  Coast  Shock  Test  Facility.  The  outputs  from  the 
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It  is  strongly  recommended  that  future  modal  testing 
be  conducted  before  the  bolt  holes  are  drilled  in  the  frame. 
An  attempt  might  also  be  made  to  weight  the  frame  with  sand¬ 
bags  to  simulate  as  nearly  as  possible  the  clamped  boundary 
conditions  input  to  the  code. 
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values  from  the  modal  testing,  the  greatest  relative  error 
(using  the  modal  test  value  as  a  base)  being  6.32 %  for  the 
tee-stiffened  plate  and  2.16%  for  the  rectangular-st if fened 
plate,  and  the  average  relative  errors  being  3.16%  and  1.12% 
for  the  tee-  and  rectangular-stif fened  plates,  respectively. 

Less  promising  is  the  seeming  failure  of  the  code  to 
predict  the  fundamental  frequency  for  either  of  the  plates. 
It  is  suggested,  however,  that  this  is  not  due  to  any  errors 
in  the  code  or  in  the  inputs  to  the  code,  but  rather  is 
attributable  to  the  fact  that  the  boundary  conditions 
actually  tested  were  not  the  clamped-clamped  conditions 
specified  in  the  STAGS  inputs.  The  H-P  5451C  Fourier 
Analyzer  appears  to  be  quite  sensitive  to  relatively 
small  changes  in  the  boundary  conditions.  When  the  tests 
were  conducted,  the  plates  were  laid  face-downward  upon  a 
heavy  marble  slab.  No  attempt  was  made  to  further  clamp 
the  heavy  aluminum  frame  into  .which  the  test  plates  are 


amount  of  CPU  time  allowed  for  eigenvalue  analysis  (NSEC  on 
the  D-2  record)  seems  to  affect  that  number.  Computation  of 
eigenvalues  over  a  range  without  any  seeding  of  known  values 
is  a  slow  process  prone  to  convergence  problems.  It  should 
therefore  only  be  used  to  derive  a  value  which  can  be  used 
for  the  second  run;  or  it  can  be  eliminated  altogether  if 
one  of  the  natural  frequencies  is  already  known.  If  such  a 
frequency  has  been  determined,  it  can  be  entered  directly 
into  the  code  and  all  frequencies  closest  to  multiples  of  a 
fixed  interval  (SHIFT  on  the  D-3  record)  will  be  computed. 
This  process  is  reasonably  fast  and  will  produce  all 
existing  frequencies  up  to  the  number  (NEIG)  specified. 
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Table  V 


Modal  Comparisons  for  Rect angular-St i f fened  Plate 

Measured  Resonant  STAGS2  Predicted 

Frequency  (Hz)  Modal  Frequecv  (Hz) 


516.5652 

804.2715 

790.0587 

924.8018 

980.5009 

981.1315 

1001.6130 

1118.0952 

1204.4885 

1229.9939 

1365.3230 

1428.2920 

1651.2498 

1651.2747 

1115 . 3547 

1852.1162 

1812 .0698 

2040.0291 

2048.0015 

2154.0908 

2368.5132 

2603.9736 

2725.8828 

2176.5366 

3039.6057 


Table  IV 

Modal  Comparisons  for  Tee-Stiffened  Plate 


Measured  Resonant 
Frequency  (Hz) 


STAGS2  Predicted 
Modal  Frequency  (Hz) 


629.7410 

653.2053 

721.1747 

772.3517 

772.7611 

793.0253 

1022.4708 

957.8445 

1310.2979 

1447.0071 

1448.4597 

1480.9890 

1499.5906 

1551.0845 

1556.0732 

1590.2947 

1747.9771 

1994.4626 

2143.0918 

2043.4048 

2236.9248 

2200.9004 

2259.6445 

2275.5000 

2430.0249 

2472.2305 

2571.7832 

2629.4729 

combined  STAGS1/STAGS2  processor  without  reference  to  the 

USA  add-on.  The  first  run  simply  specifies  a  broad  range 

within  which  the  computer  can  search  for  eigenvalues;  the 

second  uses  one  of  the  mid-range  frequencies  determined  from 

the  first  run  eigenvalues  as  a  "seed"  for  further  iterations 
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in  order  to  produce  a  refined  product.  The  actual  perform 
ance  of  these  functions  is  controlled  by  the  D-2  and  D-3 
records  ofthe  STAGS  input  deck.  Since  the  procedure  is 
outlined  clearly  in  Reference  13  and  example  input  files  are 
provided  for  both  the  tee-stiffened  and  rectangular  stiff¬ 
ened  plates  in  Appendix  C,  further  elaboration  will  not  be 
provided  here. 

Modal  testing  of  the  milled  plates  was  conducted  in  the 
manner  described  in  Section  D  of  Chapter  IV.  The  results  of 
both  that  testing  and  the  STAGS2  estimates  are  provided  for 
comparison  in  Tables  IV  and  V  tee-stiffened  and  rectangular- 
stiffened  plates,  respectively. 

It  will  be  noted  that,  for  each  of  these  plates,  the 
STAGS2  predicted  values  seem  to  correspond  very  well  with 

45 

Generally  referred  to  here  as  a  STAGS2  run. 
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Since,  in  this  case,  modal  testing  of  the  milled 
plates  was  conducted  prior  to  the  computer  runs,  essentially 
the  same  frequency  range  was  used  for  the  computer  runs  as 
for  the  testing.  Over  this  range,  the  first  STAGS2  run  will 
compute  five  modes  with  their  corresponding  natural  frequen¬ 
cies.  Although  the  user's  manual  states  that  calculation  of 
up  to  twenty  eigenvalues  can  be  specified,  and  that  all 
existing  eigenvalues  up  to  that  number  will  be  computed, 
five  seems  to  be  the  true  maximum.  No  increase  in  the 


VI.  RESULTS  AND  CONCLUSIONS 


A.  MODAL  COMPARISONS 

As  stated  in  chapter  four,  the  primary  purpose  of  modal 
testing  is  to  provide  a  check  of  the  finite  element  model. 

A  properly  constructed  model  should  exhibit  substantially 
the  same  fundamental  frequency  as  that  determined  by  modal 
testing  and  should  demonstrate  close  coincidence  with  other 
frequency  values  across  the  range  tested.  Exact  concurrence 
cannot  be  expected  since  the  computer  code  is  analyzing  an 
idealized  model  which  will  not  be  affected  by  the  externali¬ 
ties  inherent  in  the  actual  test.  Also,  the  number  of 
frequencies  determined  by  the  computer  is  a  function  of  the 
number  of  iterations^than  can  be  made  over  a  given  time 
span.  Certain  economies  in  the  use  of  computer  time  must 
therefore  be  observed.  While  the  code  might  be  able  to 
calculate  the  frequencies  for  a  large  number  of  modes;  such 
calculation  will  be  purchased  with  a  heavy  expenditure  in 
CPU  time  which  is  not  really  necessary  for  an  adequate 
check . 

Accordingly,  the  scheme  followed  here  has  been  to  limit 
the  number  of  modal  frequencies  computed  by  the  code  to 
those  that  could  be  determined  in  two  runs  for  each  plate. 
Each  of  these  runs  can  be  performed  entirely  within  the 


weight.  As  has  already  been,  stated,  geometric  similarities 
between  these  two  plates  permitted  broad  generalizations  to 
be  made  between  them  even  if  exact  comparisons  could  not  be 
drawn.  The  results  of  that  run  are  shown  in  Table  III. 

No  similar  comparisons  could  be  made  with  the  plate 
with  the  longitudinal  rectangular  stiffener.  Strain  gage 
calibrations  were  therefore  arbitrarily  set  at  a  value  at 
the  higher  end  of  the  strain  gage  range  which  could  be 
expected  to  exceed  all  requirements.  In  this  case,  that 
value  was  150,000  microstrain. 

The  results  of  these  tests  will  be  discussed  in  the 
next  chapter. 


Table  III 

Calibration  Data  for  Tee-St if fened  Plate 


Gage 

EPSA  Strain* 

1 

6210 

2 

6210 

3 

17600 

4 

8500 

5 

21700 

6 

69400 

7 

69400 

8 

69400 

9 

100000 

10 

49500 

•All  values 

in  microstrain 

Gage  EPSA  Strain* 


11 

76800 

12 

7610 

13 

7610 

14 

8510 

15 

7610 

16 

21700 

17 

6210 

18 

6210 

19 

17600 

20 

58200 
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14.03  f t .  This 


much  larger.  From  Equation  2.6,  A  = 

max 

gives  a  charge  depth  to  bubble  radius  ratio  of 
14.03  n  *  °-285 

which  implies  that  some  pressure  release  will  take  place 
into  the  atmostphere. 

For  the  second  test  shot  ,  the  problem  was  quite 
different.  As  was  discussed  in  Chapter  III,  the  tee- 
stiffened  model  appeared  too  rigid  to  successfully  demon¬ 
strate  stiffener  tripping.  The  plate  with  the  longitudinal 
rectangular  stiffener  was  a  largely  tentative  attempt  to 
produce  that  action.  It  was  desired  to  view  the  responses 
of  the  stiffener  under  a  variety  of  loadings,  beginning  with 
small  charges  and  concluding  .with  ones  in  the  8  to  9£  pound 
range.  Again  use  was  made  of  what  was  available.  In  this 
case,  the  available  charge  weight  was  £  pound. 

For  a  £  pound  charge,  the  maximum  bubble  radius  is  3.02 
feet;  which  leads  to  a  charge  depth  to  bubble  radius  ratio 
of  1.33.  This  implies  that  bubble  pulsation  will  be  a 
problem.  Fortunately,  however,  the  time  of  the  pulse  is 
relatively  late,  t  =  0.17  seconds.  Thus,  by  limiting  the 
period  of  study  to  early  times  the  bubble  pulse  can  be 
effectively  neglected. 

Calibration  levels  for  the  test  shot  on  the  te- 
stiffened  plate  were  determined  by  making  an  Et  r  .. 
plate  model  used  in  Reference  3  but  with  a  9‘  p<  an  .  .  ■  _ 
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to  the  polarity  of  the  strain  gage  connections.  T'r.  , 
histories  recorded  can  be  the  exact  inverse  of  the  true  r.c.~- 
tory.  This  can  be  especially  debilitating  when  comparison.-, 
of  the  shapes  of  strain  curves  is  of  particular  importance. 

Even  more  debilitating,  but  far  less  obvious,  was  the 
•last  of  the  problems  noted. 

The  amplifiers  on  the  Ampex  FR  1300  tape  recorders  used 
to  record  the  strain  signals  are  calibrated  for  a  40%  devia¬ 
tion  from  a  one  volt  rms  input  signal,  peak-to-peak .  The 
maximum  signal  deflection  that  is  theoretically  possible  is 
about  1.4  volts  in  either  the  positive  or  negative  direction 
The  calibration  levels  input  to  the  system  are  based  upon 
this  maximum;  signals  greater  than  1.4  volts  will  saturate 
the  amplifier  and  the  playback  will  simply  be  a  straight 
line  at  the  maxima.  Yet  on  the  recordings  for  virtually 
all  of  the  strain  gages  which  survived  to  produce  intel¬ 
ligible  records,  the  output  was  greater  than  the  1.4  volt 
limit,  but  no  concurrent  indication  of  amplifier  saturation. 
Figures  6.2  and  6.3  are  offered  as  typical  representations 
of  this  phenomena  for  each  of  the  two  plates. 

A  number  of  possible  sources  for  the  strain  gage  error 
were  considered: 

1.  The  calibration  signals  were  improperly  set. 

2.  The  signals  were  being  altered  either  as  they  were 
being  input  to  or  output  from  the  Fourier  Analyzer. 

3.  The  playback  amplifiers  on  the  tape  recorders  were 
out  of  calibration. 
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4.  The  calibration  signal  was  being  input  to  the  signal 
conditioners  at  the  wrong  location. 

The  first  three  of  these  postulated  sources  of  error  were 

readily  discarded.  The  calibration  signals  recorded  with 

the  strain  gage  signals  were  precisely  those  set  into  the 

system  prior  to  each  shot;  thus  ruling  out  the  possibility 

that  either  the  calibration  levels  were  incorrect  or  that 

the  playback  amplifiers  were  out  of  calibration.  To  check 

the  second  point,  a  series  of  calibrated  sine  waves  were 

input  to  the  Fourier  Analyzer  using  the  same  connections  and 

equipment  settings  used  for  the  strain  gage  data  recordings . 

These  calibrated  signals  were  recorded  on  disk  and  then 

output  to  the  plotter  in  exactly  the  same  manner  as  the 

strain  gage  signals  had  been.  In  no  case  did  the  Fourier 

Analyzer  alter  the  calibration  of  the  signal. 

Whether  the  remaining  suggested  cause  was,  indeed,  the 
source  of  error  or  whether  .he  actual  cause  is  something 
entirely  different  cannot  be  determined  at  this  point. 
Whatever  the  actual  cause,  the  effect  remains  the  same;  the 
strain  magnitudes  derived  from  the  two  underwater  tests 
cannot  be  used  for  comparison  with  the  results  predicted  by 
the  USA-STAGS  code. 

Although  the  magnitudes  of  the  test  derived  strain  data 
are  clearly  unuseable,  the  general  shape  of  the  strain 
curves  should  be  essentially  unaffected.  In  an  attempt  to 
resurrect  the  comparative  nature  of  this  study,  it  was 


determined  that  the  shape  of  the  curves  from  the  tests 
should  be  contrasted  with  the  strain  curves  predicted  by  th 
USA-STAGS  code.  This  attempt,  too,  was  doomed  to 
f rustrat ion . 

Figures  6.4,  6.6,  and  6.8  are  representative  samples  of 

the  strain  data  calculated  by  USA-STAGS  functioning  as  a 

post-shot  analyzer  of  the  tee-stiffened  plate.  Figures  6.5 

49 

6.7,  and  6.9  are  the  corresponding  strain  gage  outputs. 

Similarly,  Figures  6.10,  6.12,  and  6.14  are  USA-STAGS 
outputs  for  the  rectangular-st if f ened  plate;  and  Figures 
6.11,  6.13,  and  6.15  are  their  corresponding  strain  gage 
histories . 

While  similarities  exist  at  some  points  between  the 

actual  strain  histories  and  the  USA-STAGS  predicted  results 
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in  general  there  is  very  little  correlation.  This  was 
unexpected  since  earlier  discussions  with  personnel  at 
Lockheed  had  indicated  that  extremely  high  degrees  of  corre 
lation  had  been  obtained  when  the  code  results  had  been 
compared  to  with  tests  conducted  by  the  Defense  Nuclear 
Agency  and  others. 

49 

Unfortunately,  after  all  of  the  problems  encountered 
with  the  strain  gages,  there  were  not  nearly  as  many  of 
these  as  might  have  been  desired. 
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Once  again  it  should  be  emphasized  that  the  shapes 
alone  of  the  curves  can  be  compared.  As  already  explained, 
the  magnitudes  of  the  strain  gage  records  are  not  correct. 
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Figure  6.4  Tee-Stiffened  Plate,  Gage  2  (USA-STAGS) 
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Figure  6.8  Tee-Stiffened  Plate,  Gage  11  (USA-STAGS) 
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Any  number  of  possible  solutions  exist  for  this  perform¬ 
ance,  of  course;  including  the  extremes  of  either  USA-STAC-S 
or  the  strain  histories  being  totally  wrong.  In  this  case, 
however,  utilizing  the  extended  capabilities  of  the  code 
itself  may  have  provided  the  most  plausible  solution. 

The  USA  postprocessor,  POSTPR ,  has  the  capability  of 
calculating  the  fluid  pressures  at  the  fluid  structural 

interface.  These  results  may  then  be  displayed  in  either 

5  ^ 

tabular  or  graphical  form.  When  this  was  done  for  the 

cases  studied  here,  a  very  interesting  result  was  obtained. 

Figures  6.16  and  6.17  represent  the  predicted  fluid 

pressures  at  the  center  of  the  tee-stiffened  and  rectangular 

stiffened  plated,  respectively.  As  can  be  seen,  USA-STAGS 

predicts  that  intense  negat ive  pressures  will  be  generated 

within  a  very  short  time  after  the  arrival  of  the  shock 

wave.  That  is,  USA-STAGS  is  predicting  massive  local 
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cavitation  across  the  faces  of  these  plates. 

Viewed  in  that  light,  the  reason  for  the  seemingly  poor 
performance  of  USA-STAGS  in  predicting  the  strain  responses 
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As  noted  earlier,  no  graphical  capability  exists  on 
the  VAX/VMS  system  at  the  Naval  Postgraduate  School.  These 
post-shot  analysis  problems  were  run  at  Lockheed  Research, 
Palo  Alto.  It  was  they  who  most  generously  provided  the 
accompanying  plots. 
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In  actuality,  local  cavitation  would  be  marked  by  near 
zero  pressures  as  the  water  flashed  into  steam  to  fill  the 
void  created  by  the  sudden  displacement  of  the  plate.  Since 
USA-STAGS  assumes  a  homogeneous  medium  this  must  necessarily 
be  represented  by  negative  pressures. 
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THE  USA  SOURCE  TAPE 


VNOEn^ATSP  3 H‘C  1  r  -NM.  ■  -  .  S  i  LSr- 

VAx  Ta  PE  LOAD  ING  I N  C  —  ..  ■  c  T I  C’*~- 

miS  ;-u»  'C"'tjlp'5  ire  iuC*  :  •  ?!  el  'J* e S  .  ;  l-.r.jr:;-, 

?'?cut3c  Iss.  ana  manuals  ’•or  a  unot'1  o  *  ‘J  ~  -■  processors  All  flies  r  a .  e 
oeen  created  with  COPY  commanas  ana  <;in  os  -e-rievos  in  t  n «  same  manner 
with  on iy  a  Pew  lines  of  input  at  the  terminal  • ia»i ration  of  tn«  tao? 
j  irt:  ;o*4  will  reveal  t  n  a  t  tne  Piles  nave  •  a  ~  i  o  u  s  tide  n  u  m  c  e  r  s  ana  mis 
is  tne  «  e  y  to  Keeping  tnose  Piles  together  -.-.at  Belong  to  a  particular 
processor  The  following  taole  contains  tn n  l a en t i P i : a t l on  code 

USA  PROCESSOR  CYCLE 


PLUMAS  1 

AUGMAT  2 

TIM  I  NT  3 

POST=P  4 

MODL 13  5 

To  load  the  TIMINT  Piles.  Por  example.  go  r;  -ne  directory  where  meg 
are  to  reside  and  simplg  da  a  COPY  TAP'S.*  -  C  C]  command  where  TAPE  is 
the  internal  name  that  has  been  assigned  t;  the  magnetic  tape  unit  This 
process  is  simplg  repeated  until  all  aesireo  directories  have  oeen 
loaned 

The  3  character  extension'  on  the  file  names  denotes  tne  tape  of 
Pile  according  to  the  following  table. 

EXTENSION  TYPE 


.  POP 
.  INC 
.  UAX 
.  OLB 
COM 
EXE 
.  UAL 


INCLUDE  PROCEDURES 
INCLUDE  PROCEDURES 
SYM30LICS 
OBJECT  LIBRARIES 
LINK  STATEMENTS 
EXECUTAQLES 
MANUALS 


for  example.  AUGMAN  UAL  is  a  Pile  containing  the  AUGMAT  manual  while 
LINKUS*S  COM  is  the  LINK  statement  for  USA-STAGS  It  should  be  noted 
that  all  LINK  statements  must  be  modified  slightly  as  the  directories 
referenced  in  them  reflect  names  that  are  in  use  here  at  Lockheed  and 
contain  users  names  that  will  not  exist  at  other  VAX  installations 
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STAGS’.  COM 

*  .-SSI  ON  CASE!  Inp  FOR  GO 5 

*  ASSIGN  CaSEI  u'-’T  ruRCOe 

*  ASSIGN  CASE!  MOD  FGRC’21 

*  RUN  STAGS  1  SXt 

*  print  case:  out 

G  STAGS.  COfl  (Sue  jti  STaGSI  anc  S“aGS2> 

*  ASSIGN  CASEI  INP  POP 00 5 

4  ASSIGN  CASE*. OUT  FORGO© 

*  ASSIGN  CASE1  MOD  FQRG21  (optional  if  modal  plots  oesireoi 

*  RUN  STAGSI  EXE 

*  ASSIGN  CASE  1  RST  F0RG2G  (previous  solution  oata<  veqd  for  restart 

*  ASSIGN  C ASE ! . SOD  F0RG22  (optional  if  solution  data  snouia  os  -svsd 

*  RUN  STAG52  EXE 

*  PRINT  CASE!  OUT 

N.  POST?  COM  (Eiecute  STAGSI  and  PQS”P ) 

*  ASSIGN  CASE*  INP  F0RGC5 

*  ASSIGN  CASE*  OUT  FOP 006 

*  RUN  STAGSI  EXE 

*  ASSIGN  CASE 4  PPN  F0R0G5 

*  ASSIGN  CASE4  POP  F0RG06 

*  ASSIGN  CA5E4  RST  F0R020 

*  ASSIGN  CASE 4  SOD  FGRC22 

*  RUN  POSTP  EXE 

*  DEL  FGR009. DAT 

I.  STAPL.COM  ( E  i  acute  STAPD 

*  ASSIGN  CASE l  PIN  FOR 005 

*  ASSIGN  CASE  1  RST  F0PC20 

*  ASSIGN  CASE l . MOD  F0R021 

*  ASS I CN  CASE 1  OUT  FGRG06 

*  RUN  STAPL. EXE 

*  (Eiecute  plot  file  on  local  plot  davica! 


(May  ba  omitted  if  F0R002  DAT  saved/ 
(previous  solution  data  from  CASE4  SOD) 
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PROCEDURES 


A.  L13R  COM 


;  i  «  S'-ograms  me  c  r 


l  ~QP  STAGSU  foR/nOlIST 
t  l I  UP AR  y /CREATE  STaGSu  OLB  STAGSU 

*  MACRO  8LI0MF  MAS 

*  MACRO  J03STA  MAP 

t  LIB  STAG 30  BUOMr,  UG3STA 

»  FOP  STAG31.  FOR/NOLIST  _ 

*  LIBFAPY/CREATE  STAGS1  0L3  5TAGS1 

*  FOR  MQUNTEL.  FOR/NOLIST 

*  LIB  STAGS l  MOUNTEL 

*  FOP  3TAGS2  FOR/NGlIST 

*  LIBRARY/CREATE  5TAGS2  OLB  STAGS2 

*  FOR  MOUNTEL.  FOR /NOLI ST 

*  LIB  STACS2  MOUNTEL 

*  FOR  STAPL  FOR/NQLIST 

*  i_IBPARY/CR£AT£  STAPL  OLB  STaPL 

*  POP  POSTP  FOR / NOLIST 

*  LIBRARY/CPEATE  POSTP  OLB  POSTP 

*  FOR  TP  FOR /NOLI ST 

*  FOR  ZSYS  F0P/N0LI3T 

*  LIB/CREATE  ZSYS  ZSYS 

*  FOR  RES. FOR/ NOL 1ST 

B.  LSTACS1.COM  (Load  STAGS  1 > 

*  LINK /MAP -STAGS l /BP IEF/ EXECUTABLE-STAGS  1- 

STAGSl  OLB/LIB/ INCLU0E=MAIN1, - 
STACSU.  OLB/LIB 

C.  LSTAGS2.COM  (Load  STAGS2 ) 

*  LINK/MAP-STAGS2/BRIEF/EXECUTABLE-STAGS2- 

STAC32.  0LS/LIB/INCLUDE»MAIN2.  - 
STAGSU  OLB/LIB 

D.  LSTAPL.COM  (Load  STAPL) 

*  LINK/MAP-STAPL/8RIEF/EXECUTA0LE-STAPL- 

STAPL.  OLB/LIB/ INCLUDE-3TAPL.  - 
STAGSU  0L8/LIB.  - 

(Fil*  containing  local  plot  routine*) 

E.  LPOSTP  COM  (Load  POSTP) 


»  LINK /MAP -POSTP  /  BRIEF  '  E/EC'Jt  ABLE=POSTP- 
POSTP  OLB/LIB/ InClUOE-PGSTP,  - 
STAGS1  OL8/LI3.  STAGS2  OLB/LIB. - 
3TAGSU  GLS/LI3.  NICE.  OLB/LIB 


A.  THE  STAGS84  SOURCE  TAPE 


-£r.£'Cli  Au.,037  i9s4 


5TAGSC  l  PS  OSS  am  Tape 
VAX /VMS  OPERATING  Sf'5TE" 


1  TAPE  FORMAT 

9-T  r«ck 
1600  BP  I 

Labeled  (STAG34) 


2  TAPE  CONTENTS 


lie 

Name 

Tgpe 

Description 

1 

INFO  VAX 

Te « t 

Table  of 

Cont?r,rj 

*  proc eiurf? 

2 

STACSU  FOR 

Source 

FORTRAN 

prog 

for  STAGS  utility  library 

3 

3LICMF  MAR 

Source 

MACRO  program  *cr 

STAGS  utility  library 

4 

J0G3TA  MAR 

S  o  u  r  c  r 

MACRO  program  ?or 

STAGS  utility  library. 

5 

STAGS  1  FOR 

Source 

FOR  THAN 

prog  * iin i 

forSTAGSl  library 

6 

MOUNTEL. FOR 

Source 

FORTRAN 

progs*  am« 

Mr  th*  MOUNT  elements 

6 

STAG 32  FOR 

Sourc  e 

FORTRAN 

programs 

for  5TAGS2  library 

S 

STAPL.  for 

Source 

FORTRAN 

prog  - 

for  $tAPL  library 

9 

POSTS  FOR 

Source 

FORTRAN 

p  r  og ” am* 

for  POSTF  library 

10 

STAGS  1  OLB 

Liorarg 

L  i  br arg 

or  Sr-C-rl 

programs 

11 

STACS2  OLB 

L  l  b  r  ar  g 

L  l  b  r  a  r  g 

of  ST-GS2 

program*!. 

12 

STAG3U  OLB 

Li  brarg 

Librarg 

or  STAGS 

utility  routines. 

13 

POSTP  OLB 

Li  b  rar  g 

Liorarg 

of  P-jZT? 

programs 

14 

STAPL.  OLB 

L  l  b  r  ar  g 

L  i  brarg 

of  STAPL 

plot  prog  rams. 

15 

NICE  OLS 

L  i  brarg 

16 

TP  FOR 

Source 

Thurston 

Processor. 

17 

TP.  OBJ 

Object 

Thurston 

Proce  *  s  or . 

IS 

ZSYS.  FOR 

Source 

Z-sgstem 

i  utilities  (for  TP  > . 

19 

RES.  FOR 

Source 

Reduced 

Equa  t i on 

Solver 

20 

RES  OBJ 

Object 

Reduced 

Equa  1 1 on 

Solver 

21 

ZSYS  OLB 

Librarg 

Z-sgstem 

utilities  object  librarg 

!-36 

CASE*.  *.  * 

List 

Input  and  output 

for  5  samp  le  cases. 

3.  LOADING  FROM  TAPE 

Use  the  following  control  cord*  to  copg  oil  files  to  due 

*  MOUNT  MTAO:  STAGB4  TAPE 

*  COPY  TAPE:  *.  *.  *  n 

File*  mag  also  b»  read  indiduallg  from  tape  with  the 
fallowing  control  card*  <e  j.  .  to  load  STAGS!  OLD i . 

*  COPY  TAPE: STAGS  1. 0L8  C] 

Complete  sample  case  output  mag  be  oorained  as  follows 

*  PP INT / HEAD  CASE*  *  * 
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APPENDIX  A 


SOURCE  DATA 

This  appendix  contains  the  information  listings  supplied 
with  the  STAGS  and  USA  source  tapes.  Section  A  is  a  full 
listing  of  the  contents  of  the  STAGSS4  tape  along  with 
instructions  for  mounting  the  tape,  command  files  for 
creating  and  linking  executables,  and  sample  command  files 
for  running  the  code.  Section  B  contains  similar  listings 
for  the  USA  source  tape  excepting  that  no  command  files  are 
listed.  USA  places  the  link  commands  on  the  tape  itself  for 
immediate  access.  Command  files  for  running  the  modules  are 
identical  in  construction  to  those  for  STAGS  and  can  be 
adapted  from  those  in  Section  A. 


In  particular  the  question  of  whether  the  local  cavita¬ 
tion  predicted  by  USA-STAGS  is  actually  occurring  or  is 
simply  a  result  of  the  choice  of  model  must  be  addressed. 

It  is  recommended  that  this  study  should  be  redone,  but 
with  the  backing  structure  included  as  a  part  of  the  basic 
model.  POSTPR  could  then  be  run  to  determine  whether  cavi¬ 
tation  was  still  occurring  across  the  plate.  If  cavitation 
was  still  being  indicated,  then  either  a  new  model  would 
have  to  be  found  or  the  cavitated  region  would  have  to  be 
modelled  using  a  code  capable  of  handling  it. 

Whatever  the  findings  of  continued  tests,  it  may  have  to 
be  recognized  that  increasing  the  scope  of  the  models 
studied  may  put  USA-STAGS  beyond  the  reach  of  masters-level 
research  at  the  Naval  Postgraduate  School.  USA-STAGS  is  a 
diverse,  highly  capable  code  which  has  the  reputation  for 
high  levels  of  accuracy.  The  price  that  is  paid  for  such 
capability,  however,  is  size  and  complexity.  Given  the 
limited  time  and  computer  resources  available  to  most 
students  at  NPS ,  USA-STAGS  may  be  beyond  their  capacity.  If 
such  research  is  undertaken,  it  is  most  emphatically  recom¬ 
mended  that  a  1 '.rge  block  of  time,  weeks,  perhaps,  be  spent 
at  Lockheed  Research  in  Palo  Alto,  California  learning  the 
code  and  perfecting  models. 


1.  Rather  than  using  calibration  resistors,  the  strain 
gage  bridge  should  be  balanced  using  strain  gages 
identical  to  those  on  the  test  plate.  These  should 
be  attached  to  an  aluminum  plate  in  the  terminal  box. 
This  has  two  effects:  First,  the  balancing  gages 
would  be  subjected  to  the  same  changes  in  temperature 
as  the  plate  gages.  Second,  the  balancing  leg  of 

the  bridge  would  be  of  the  same  length  as  the  measure¬ 
ment  leg.  Combined,  these  would  make  the  bridge 
far  less  sensitive  to  environmental  changes. '  A 
zero  balance  could  be  maintained  and  calibration 
fluctuations  would  be  less  likely. 

2.  A  zero-time  triggering  signal  should  be  superimposed 
upon  all  instrument  channels..  This  would  permit  an 
exact  assessment  of  the  time  of  initial  plate  response 
to  the  incident  wave. 

3.  Further  shielding  and  grounding  of  the  pressure  gage 
cables  should  be  attempted  to  prevent  the  saturation  of 
the  pressure  gage  amplifiers  through  currents  induced 
by  charge  detonation. 

From  the  results  that  have  been  explored  in  the 
preceding  section  of  this  chapter  it  is  clear  that  certain 
questions  need  to  be  raised  about  the  validity  of  the  flat 
plate  model  that  has  been  used  here.  These  same  questions 
should  be  asked  about  the  ability  of  USA-STAGS  or  any  other 
code  to  handle  that  model  under  the  conditions  imposed. 


141 


the  eight  pound  TNT  charge  as  the  standard  ■  [Ref .  3]  is  still 
valid.  It  is  therefore  recommended  that  that  weight  charge 
be  maintained  if  at  all  possible. 

Second,  the  lack  of  adequate  strain  gage  data  has  been 
an  incapacitating  flaw  in  the  test  results.  Whether  the 
poor  quality  of  that  data  has  been  the  result  of  errors  in 
the  equipment  set-up  at  the  West  Coast  Shock  Test  Facility, 
the  interface  between  that  equipment  and  the  strain  gages  on 
the  test  plate,  or  some  other  cause  is  uncertain.  What  .is 
certain  is  that  the  present  arrangement  of  sending  the  test 
plate  and  backing  structure  to  WCSF  a  few  days  before  the 
test  shot  for  attachment  and  then  allowing  only  a  few  hours 
immediately  before  the  shot  for  equipment  testing  is 
insufficient.  Development  of  an  in-house  instrumentation 
capability  at  the  Naval  Postgraduate  School  is  considered  to 
be  essential.  This  would  permit  the  gages  and  electronic 
equipment  to  be  integrated  into  a  single  system.  Thorough 
testing  could  be  conducted  before  the  actual  date  of  the 
test;  and,  equally  important,  responsibility  for  the  outcome 
of  the  test  would  reside  in  a  single  organization. 

Changes  in  the  electronic  instrumentation  have  already 
been  recommended  [Ref.  3].  Enactment  of  those  changes  would 
undoubtedly  improve  the  data  responses.  A  few  other  modifi¬ 
cations  suggested  by  the  results  of  this  study  are  consid¬ 
ered  equally  important,  however. 
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the  second  objective  presented  in  the  introduction,  that  of 
providing  guidelines  for  future  study,  that  definite  contri¬ 
butions  can  be  made. 

In  the  realm  of  underwater  explosion  testing,  some 
specific  problem  areas  were  noted  in  the  preceding  chapters. 
Solutions  to  some  of  these  (e.g.,  strain  gage  choice  and 
attachment)  were  presented  in  the  text  of  those  chapters. 
Others  of  these  problems  are  of  a  more  fundamental  nature 
and  must  be  discussed  separately. 

First,  a  satisfactory  plan  of  testing  cannot  be  devel¬ 
oped  without  an  adequate  source  of  appropriately  sized 
explosive  charges.  The  inability  to  match  charge  weight  to 
test  platform  dimensions  makes  it  impossible  to  use  any 
computer  code  in  its  preshot  capacity  and  impacts  negatively 
upon  the  reproduceability  of  test  results.  Either  the 
charge  weight  must  be  standardized  or  the  test  platform 
dimensions  must  be  altered  to  accommodate  changes  in  charge 
weight.  In  either  case,  the  most  apparent  solution  seems  to 
be  a  carefully  orchestrated  test  plan  with  long  lead  times 
before  each  test  shot.  If  a  match  is  not  possible,  it  is 
suggested  that  the  shot  should  probably  be  aborted.  That 
valuable  information  was  obtained  from  such  mismatched 
conditions  in  this  study  is  considered  to  be  more  a  charac¬ 
teristic  of  a  program  in  its  developmental  stages  than  of 
well-conceived  experimental  practice.  Until  more  experience 
is  gained,  it  appears  that  the  rationale  dictating  use  of 
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force  of  the  incident  pressure  wave  may  not  ,  in  fact  ,  'oe 
impinging  upon  the  test  plate;  since,  while  that  wave  is 
"pushing"  on  the  front  side  of  the  structure,  it  is  also 
"pulling"  on  the  back. 

Whether  this  decrease  in  force  would  be  sufficient  to 
eliminate  cavitation  is  an  open  question.  More  testing  must 
be  done.  For  the  USA-STAGS  code  to  be  able  to  predict  this, 
however,  it  is  not  enough  to  model  only  the  test  panel.  The 
entire  assembly  of  test  panel,  milled  edge,  and  backing 
structure  must  be  included.  If  local  cavitation  is  still 
predicted,  then  additional  action  will  be  required  because 
Doubly  Asymptotic  Approximation  which  is  the  basis  for 
the  USA  code  is  not  valid  where  discontinuities  such  as 
cavitation  exist  between  the  fluid  and  the  structural 
surface55  [Ref.  17]. 

C.  CONCLUSIONS  AND  SUGGESTIONS  FOR  FUTURE  STUDY 

Clearly,  problems  encountered  both  in  conducting  under¬ 
water  testing  and  in  using  the  USA-STAGS  code  have  precluded 
any  firm  conclusions  from  being  made  about  the  predictive 
capabilities  of  the  code.  Thus,  the  primary  objective  of 
this  study  has  not  been  met.  It  is  in  the  fulfillment  of 

55 

Lockheed  has  developed  a  processor,  USA-STAr '  TA 
(Cavitated  Fluid  Analyzer)  which  can  be  used  to  model  the 
cavitated  region.  The  process  that  must  be  followed  is  to 
model  layers  of  fluid  at  increasing  intervals  from  the  plate 
until  no  further  cavitation  has  been  found.  At  that  point, 
the  standard  DAA  approximation  may  be  applied. 
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discontinuity  is  therefore  being  created  in  the  fluid 
structural  interface.  This  discontinuity  is  manifested  by 
the  pressure  drop  which  causes  the  local  cavitation. 

The  strain  gage  histories  seem  to  imply  that  local  cavi¬ 
tation  is  indeed  occurring.  Whether  that  is  in  fact  the 
case,  or  whether  the  cavitation  predicted  by  USA-STAGS  is 
simply  a  function  of  the  way  that  the  plates  are  modelled  is 
something  that  must  still  be  determined,  however. 

In  these  studies,  the  plate  alone  has  been  modelled. 

The  structure  that  is  ''known"  to  the  computer  code  can 
therefore  be  conceived  as  infinite  in  extent  and  must  neces¬ 
sarily  bear  the  full  brunt  of  the  pressure  wave.  In  actu¬ 
ality,  this  is  not  so.  The  test  plate  is  finite  in  size  and 
pressure  relief  will  occur  around  its  boundaries.  Since  the 
breadth  of  the  shock  front  is  on  the  order  of  five  feet  for 
charges  of  the  weight  used  in  these  tests  (this  is  a  func¬ 
tion  of  the  speed  of  sound  in  the  fluid  and  the  pressure 
decay,  it  will  vary  with  charge  weight),  the  entire  assembly 
will  be  enveloped  by  the  pressure  front.  This  means  that, 
almost  instantaneously  after  the  incident  pressure  is  felt 
by  the  test  plate,  the  pressure  will  begin  to  relieve  around 
the  edge  of  the  plate  and  a  negat ive  pressure  will  be 

experienced  on  the  backs  of  the  flanges,  and,  soon  after,  on 

54 

the  rear  of  the  backing  structure.  Therefore  the  entire 
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Complex  interference  patterns  can  also  be  expected. 


can  be  better  understood,  and  even  the  strain  gage  responses 
themselves  become  more  understandable. 

For  example,  in  Figure  6.11  an  initial  compressive  strai 
can  be  observed  as  the  plate  reacts  to  the  incident  pres¬ 
sure.  This  does  not  occur  immediately  but  is  delayed  some¬ 
what  as  the  inertia  of  the  plate  is  overcome.  In  this 
instance,  the  initial  compressive  strain  is  mimicked  almost 
exactly  by  the  response  predicted  by  USA-STAGS.  The  strains 
then  return  to  zero  or  become  tensile  as  cavitation  occurs. 
As  would  be  expected,  the  responses  predicted  by  the  code 
are  much  more  violent  then  those  seen  in  the  strain  gage 
histories.  This  is  so  because  they  are  activated  by  large 
values  of  assumed  pressure  rather  than  the  actual  zero  pres¬ 
sure  that  would  occur  with  cavitation. 

The  other  plots  seem  to  confirm  this  basic  trend: 

Initial  zero  or  compressive  strain  depending  upon  the  loca¬ 
tion  of  the  element /st rain  gage  followed  by  a  positive  trend 

53 

in  the  strain  history  as  cavitation  develops. 

The  actual  mechanisms  involved  in  the  results  that  seem 
to  be  indicated  here  are  not  fully  understood .  The  implica¬ 
tion,  however,  is  that  the  velocity  of  the  plate  as  it  is 
being  deformed  is  much  greater  than  that  of  the  fluid.  A 

^It  is  suggested  that  Figures  6.9  and  6.15  also  adhere 
to  this  basic  pattern,  but  are  subject  to  the  polarity 
problem  mentioned  earlier.  It  is  believed  that  these  plots 
should  actually  be  inverted. 
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APPENDIX  B 


BASIC  SHELL  TYPES  FOR  USE  WITH  STAGS 

This  appendix  contains  a  detailed  listing  of  the  eleven 
basic  shell  types  used  with  the  STAGS  code.  Drawings  are 
included  to  show  numbering  conventions  for  edges  and 
corners.  This  information  can  be  particularly  useful  for 
discretizing  the  shell  unit  and  specifying  compatibilities 
with  .adjoining  shells.  (All  figures  from  Reference  13.) 
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PROP(')  *  x, 
PRCP(2)  =  X . 
PROP'  3  f  =  r? 
PROp(4j  s  y. 


Subscripts  refe 
corner  number. 


QUADRILATERAL  PLATE: 

*1  =  (1-XJ  •  (1-Y),  £2  sX(l.Y),  fj  s  XT. 


x(XY)  =  C  {.  X. 

i=l  1  i 

Y  (X  Y)  S  z  f.  Y. 

i-l  1  1 

PROPIU  =  X, 
PRCP(2>  r  Y 
PROP!  3)  *  X, 
prop(4)  =  y; 
PROP!  5)  =  XT 
PROP!  6)  =  y; 
PROP! 7)  r  x* 
PROP(3)  :  y^ 


ANNLXAR  PLATE. 


y  s-  z*  *  r“  s  X* 
x  *  f(X,  Y)  «  0 
y  *  g  (X,  Y)  =  X  sia  Y 
J  *  h  (X.  Y)  *  X  cos 


PROP!  i )  =  r 
PROP(’)  *  Ru 
PRCPO)  =  a 
PROP(4)  =  ^ 


Lateral  displacerwnt  posit in? 
( negative  x } ! 


APPENDIX  C 


USA-STAGS  INPUT  FILES 


This  appendix  contain? 
complete  USA-STAGS  runs  f 
rectangular -st  i  f  '  » * 

in  which  they  w  •  u ;  :  •  ~  , 

-  ST  AGS  2  \!f  it A  -  t  ■ 

-  STAGS  1  I  r.y  .■  ■ 

-  Fluid  Mas-  I  • 

-  Augmen*  er  ■' v ' 

-  T  ime  1  .n  *  -*g  -  <  * 

-  USA  p- 

The  ir.pur  ::1 «-  * 

examp  1 e  o  f  5TA1 1  : 

stiffened  p  1  a  t  “  -  '  i  .  ■ 

Since  these  las*  f  ;  :  i  •  • 

considered  necessar-.  *  n  * 


he  input  files  used  for  th 
•  *  r.  *  •>•  t  -e-s*  l  f  fened  and 

T~~  ■»-  i  *  •*']  .  r  *  he 


t  ■ 

•  ••  *  •  :gu  .  a  : 

i  :  .  .  :  *  vi.-, 

:  :*•  ♦» x amp  .  •  -  t  ■  ■  -  t  r.  .ms s- 
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TEESTAG32. IN P;  1  0 


9-SE3-19«4  1-3:5? 


STAGS?  OATA  FDR  TRANSVERSE  TEE  STIFFENER  (MODAL  ANALYSIS) 

20000010053-1 

50045  3-2 

:  0  2  $  3-3 

0.  S  C-l 

1  0  10000  S  0-2 

12  200.  0.  2600.  S  0-3 

1  1  6  1  1  1  1  1 1  3  t 1  3  1 1  3  5  F-l 

1  2  2  4  J  G-l-1  COM? A  TA8ILITY  1-2 

123«5  G-l-2  C0MPATA3ILITY  1-3 

3244*  G-l-3  COMPATASILITY  3-4 

3  2  5  2  5  G-t-4  CO-MPAT ABILITY  3-5 
1051-1 

10. Eb  .3  0.  .093  S  1-2  MATERIAL  PROPERTIES 
1115  K  —  l — 1  FLOOR  SHELLS 
1  .1875  5  K-2-l 

2115  K-l-2  STIFFENER  SHELLS 

1  .125  $  5-2-2 

2  0  $  m-i -i 

0.  fa.  0.  3.  5  M-2A-1 
100. ■ -.0938  5  M-5-t 
410  $  N-l-1 

46245  P-1-1  BOUNDS  1 

0  5  0-1-1 

0  $  R-l-1 

205  M-l-2 

0.  fa.  3.  9.  $  M-2A-? 

1  0  0.  -.0938  5  M-5-2 
410  5  N  -  1  -  2 

4  2  2  6  5  P-1-2  80UN0S  2 

0  5  0-1-2 

0  5  R-l-2 

225  M-i-3 

0.  fa.  0.  .7313  $  M-2A-3 

1  3  0.  3.  .7813  5  M-3-3 

2  0  90.  0.  5  m-5-3 
410  5  N-l-3 

4  fa  2  fa  5  P-1-3  BOUNDS  3 

0  I  0-1-3 

0  5  R-l-3 

225  M— i-4 

0.  fa.  3.  3.3594  s  M-2A-4 

1  3  0.  3.3594  .7813  5  M-3-4 
20  0.  .0625  5  m-5-u 

410  3  N-t-4 

4  3  2  fa  5  P-1-4  90'JNOS  4 
0  5  0-1-4 
0  5  R-l-u 

2  2  5  M-I -5 

0.  fa.  2 . fa40o  5.  5  m - 2  *  - 5 

1  1  0.  2.6406  .7813  I  «-3-5 

2  0  0,  .0625  3  m • 5 - 5 
410  5  N- l -5 

4  fa  2  3  5  P-1-5  90UNPS  5 
0  5  0-1-8 
0  5  R-l-5 


STIFFENER  (PLASTIC) 


STagsi  data  pop  transferee  tee 
6  1  1  t  0  0  1  0  1  *  3-1 
5004*  3-2 
102*3-3 
0  *  C-l 

0  1000  1  1  *  D-l 
0  2.  £-3  20  £-6  *  £-1 

0  4  *  £-2 

11  6  11  11  11  3  11  3  11  3  *  F-l 
l  2  2  4  *  G-l-l  COMPATAB ILITY  1-2 

1  2  3  4  *  G-l-2  COMPATAB ILITY  1-3 

3  2  4  4  *  0-1-3  C0MPATA3 ILITY  3-4 

3  2  5  2  *  S-i-4  C0MPATA3 ILITY  3-5 

131*1-1 

10. £6  .3  0.  .098  *  1-2  MATERIAL  PROPERTIES 

.0041  41.  £3  06  45. E3  1.  47. £3  s  1-3  PLASTICITY 
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